Global phase diagram of two-dimensional Dirac fermions in random potentials 
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Anderson localization is studied for two flavors of massless Dirac fermions in two-dimensional space 
perturbed by static disorder that is invariant under a chiral symmetry (chS) and a time-reversal 
symmetry (TRS) operation which, when squared, is equal either to plus or minus the identity. The 
former TRS (symmetry class BDI) can for example be realized when the Dirac fermions emerge 
from spinless fermions hopping on a two-dimensional lattice with a linear energy dispersion such as 
the honeycomb lattice (graphene) or the square lattice with 7r-flux per plaquette. The latter TRS 
is realized by the surface states of three-dimensional Z 2 -topological band insulators in symmetry 
class CII. In the phase diagram parametrized by the disorder strengths, there is an infrared stable 
line of critical points for both symmetry classes BDI and CII. Here we discuss a "global phase 
diagram" in which disordered Dirac fermion systems in all three chiral symmetry classes, AIII, CII, 
and BDI, occur in 4 quadrants, sharing one corner which represents the clean Dirac fermion limit. 
This phase diagram also includes symmetry classes All [e.g., appearing at the surface of a disordered 
three-dimensional Z 2 -topological band insulator in the spin-orbit (symplectic) symmetry class] and 
D (e.g., the random bond Ising model in two dimensions) as boundaries separating regions of the 
phase diagram belonging to the three chS classes AIII, BDI, and CII. Moreover, we argue that 
physics of Anderson localization in the CII phase can be presented in terms of a non-linear-sigma 
model (NLctM) with a Z 2 -topological term. We thereby complete the derivation of topological or 
Wess-Zumino-Novikov-Witten terms in the NLcrM description of disordered fermionic models in all 
10 symmetry classes relevant to Anderson localization in two spatial dimensions. 
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I. INTRODUCTION 

A. Dirac fermions in condensed matter physics 

Massless Dirac fermions emerge quite naturally from 
non-interacting and bipartite tight-binding Hamiltonians 
at low energies and long wave-lengths when the fermion 
spectrum of energy eigenvalues is symmetric about the 
band center and the Fermi surface reduces to a finite 
number of discrete Fermi points at the band center. This 
situation is generic for non-interacting electrons hopping 
with a uniform nearest-neighbor amplitude t along a one- 
dimensional chain. For non-interacting electrons hopping 
on higher dimensional lattices, this situation is the excep- 
tion rather than the rule, for it is only fulfilled when the 
hopping amplitudes are fine-tuned to the lattice. 

In the case of graphene, when described by the uni- 
form hopping amplitude t between the nearest-neighbor 
sites of the honeycomb lattice, there are two bands in the 
Brillouin zone of the underlying triangular Bravais lat- 
tice that touch at the 6 corners of the Brillouin zone [see 
Fig. ^(a)]. El Because the unit cell contains two sites and 
because the number of inequivalcnt Fermi points is two, 
these Dirac fermions realize a four-dimensional represen- 
tation of the Dirac equation in two-dimensional space if 
we ignore the spin degrees of freedom. 

For non-interacting spinless electrons hopping on the 
square and (hyper-)cubic lattices, Dirac fermions emerge 
in the vicinity of the band center whenever the trans- 



lation invariance of the lattice is broken by choosing 
the sign of the nearest-neighbor hopping amplitudes of 
uniform magnitude t in such a way that their prod- 
ucts along any elementary closed path (a plaquette) is 
— t 4 [see Fig. 0(b)]. This pattern of nearest- neighbor 
hopping amplitudes preserves time-reversal symmetry. 
It amounts to threading each plaquette by a magnetic 
flux of 7T or, equivalcntly, — ir in appropriate units and 
is thus called the 7r flux phase. In the 7r-flux phase 
for the d-dimcnsional hypcrcubic lattice, there are 2 d 
non-equivalent sublatticcs. Correspondingly there are 2 d 
Fermi points and the emerging Dirac Hamiltonian in the 
vicinity of these Fermi points is 2 d dimensional. Because 
the minimal irreducible representation of the Dirac equa- 
tion in d dimensions is 2^ d+1 )/ 2 l dimensional ([x] denotes 
the largest integer smaller than or equal to x), the tt- 
flux phase yields a representation of the Dirac equation 
larger than the minimal one in all dimensions except for 
d = 1 . This is called the fermion doubling problem, for it 
prevents a lattice regularization of the standard model of 
Elementary Particle Physks that represents its particle 
content (quarks, leptons). □ 

The fact that the fermion-doubling problem affects 
both graphene and the 7r-flux phase in two dimensions 
is not a coincidence. The fermion-doubling problem is 
a generic property of non-interacting local light-binding 
Hamiltonians with time-reversal symmetry!] 

It is possible to circumvent the fermion-doubling prob- 
lem in the following way. 

We consider first a one-dimensional chain along which 
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FIG. 1. (a) Hexagonal Brillouin zone of graphene with the 
conduction and valence bands touching at the zone corners 
in the linear approximation. There are 6/3 = 2 inequivalent 
Fermi points (Dirac cones), (b) The 7r-flux phase for the cubic 
lattice assigns the nearest-neighbor hopping amplitudes +t 
for the thin bonds and —t for the thick bonds with t a real 
number. There are 8 = 2 3 inequivalent sublattices labeled 1 
to 8. 



a spinless electron hops with the uniform nearest- 
neighbor amplitude t. We also impose periodic boundary 
conditions [see Fig. |](a)]. We fold the spinless electron's 
dispersion on half of its Brillouin zone and open a gap at 
the folded zone boundaries by dimerization of the hop- 
ping amplitude, t — > t ± dt, as it occurs for example 
through its interaction with an optical phonon within a 
Born-Oppenheimer approximation. At low energies, the 
effective fcrmionic Hamiltonian is the one-dimensional 
massive Dirac equation with the mass set by the dimen- 
sionless parameter St/t assumed to be smaller than unity. 
Imagine now that the dimerization pattern is defective at 
two sites that are far apart relative to the characteristic 
length scale (t/5t)a where a is the lattice spacing [see 
Fig. ||(c)]. At the level of the effective Dirac equation, 
this means that the mass term changes sign twice, once at 
each defective site. Two bound (i.e., normalizablc) states 
appear in the spectrum [see Fig. ||(d)] with the remark- 
able property that they have opposite helicity (chirality) 
and an exponentially small overlap or, equivalently, en- 
ergy splitting, for they are exponentially localized with 
the localization length .of, order (t/6t)a around their re- 
spective defective sites. QU 

The same mechanism applies in any d-dimensional 
space, be it for the massive Dirac equation,!] or for 
tight-binding Hamiltonians with sublattice symmetry 
[see Fig. H(e)],Llia an d has been used in lattice gauge 
theory as pa means to overcome the fermion doubling 
problcmoE] For example, the massive Dirac equation 
in odd d-dimcnsional space supports massless bound- 
ary states with a common helicity (chirality) along each 
even (d — l)-dimensional boundary where the mass term 
vanishes. A complete classification of all such two- 
dimensional boundary states was part of the classifi- 
cation of topological insulators in spatial dimensions 
d = 1,2,3 given in Ref. |ll| in terms of the generic sym- 



metry classes arising from the antiunitary operations of 
time reversal and particle-hole symmetry [underlying the 
work of Altland and Zirnbauer on random matrix theory 
(RMT)].E3c3 A systematic regularity (periodicity) of the 
classification as the dimensionality is varied, in general 
dimension, was discovered upon the use of K-Theory by 
Kitaevlia (see also Ref. [l6]). As shown in Refs. [Tt] and 
H this can, alternatively, be understood in terms of the 
lack of Anderson localization at the boundaries. More 
recently, an understanding of this classification of topo- 
logical insulators in terms of quantum anomalies was de- 
veloped.!!^ 



B. Anderson localization for Dirac fermions in two 
dimensions 

Anderson localization^ for non-interacting two- 
dimensional Dirac fermions was first studied in narrow 
gap semiconductors by Fradkin in 1986.EH This work 
was followed up in the 90's with non-perturbativc re- 
sults motivated by the physics of the integer quantum 
Hall effect (IQHE), the randorp^bojftd Jaing model, and 
dirty <i-wavc superconductors£j€3oT2d With the re- 
cently available transport measurements in mesoscopic 
samples of graphene, as well as the identifications 
alloy Bi 1 _ 2 ,Sb 2 . in a certain-xange of compositions x,i 
the compounds Bi 2 Te 3 pEl Sb 2 Te 3 £3 and Bi 2 Se 3 , 
and the prediction for another 50 and counting mate- 
rials as three-dimensional Z 2 -topologicaA band insulators 
that support surface Dirac fermions the localization 
properties of random Dirac fermions have become rele- 
vant from an experimental point of view. 

While all these examples share the massless Dirac spec- 
trum as the energy dispersion in the non-interacting and 
clean limit, the effects induced by randomness - weak 
localization, universal conductance fluctuations, localiza- 
tion, metal-insulator transition, spectral singularities, etc 
- vary with (i) the intrinsic symmetries respected by the 
disorder, (ii) the dimensionality of the Dirac matrices 
representing the Dirac Hamiltonian, and (iii) the strength 
and/or correlations in space of the disorder. 

When space is effectively zero-dimensional, i.e., at the 
level of RMT, ten symmetry classes have originally been 
identified and labeled according to the £a*tan classifica- 
tion of symmetric spaces (see Table ||) .E^Hlil 

As emphasized in Refs. [n]and [l2], the two-dimensional 
fcrmionic replicated NLerMs in eight of the ten sym- 
metry classes allow for terms of topological origin, in 




the form of cither 9 term; 
Wittcn (WZNW) termsM 



or Wess-Zumino-Novikov- 
(see Table |). Symmetry. 



classes A, C, and D support Pruisken (8) terms.cZrH 
Symmetry classes AIII, Dili, and CI support WZNW 
terms. Finally, symmetry classes All and CII support 
Z 2 -topological terms. 

WZNW terms in symmetry classes AIII, Dili, and CI 
appear when Dirac fermions propagate in. the presence of 
static vector-gauge-like randomness H3t21I This can only 
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TABLE I. Table of topological terms that can be added to thenreplicated fermionic non-linear-sigma model (NLcrM) describing 
Anderson localization in two dimensions and the classificationtJ of topological insulators (superconductors) in three dimensions. 
Symmetry classes indicated by the "Cartan label" are classified according to the presence or absence of time-reversal, particle- 
hole, and "sublattice" symmetries which we abbreviate as TRS, PHS, and SLS, respectively. The presence of TRS and PHS 
is denoted by or "—1," depending on whether the square of the (antiunitary) operator implementing the symmetries 

equals +1 (identity) or — 1, whereas the presence of SLS is denoted by "1". The absence of these symmetries is denoted by 
"0". The SLS is a product of TRS and PHS. For historical reasons, the first three rows of the table are also referred to as 
the orthogonal, unitary, and symplectic symmetry classes. When the disorder respects a sublattice symmetry as in the next 
three rows, the terminology chiral is also used. Finally, the last four rows can be realized as random Bogoliubov-de-Gennes 
(BdG) Hamiltonians. Target spaces for fermionic replicated NLcrM (N is the replica index and the limit N — > is understood) 
are given in the fifth column. The penultimate column lists the nature of the topological term compatible with the target 
and two-dimensional base spaces. The symbols Z and Z 2 in the last column indicate that the topologically distinct phases 
within a given symmetry class of topological insulators or superconductors in three spatial dimensions are characterized by an 
integer topological invariant (Z) or a Z 2 quantity. The symbol "0" denotes the case when there exists no topological insulator 
(superconductor) . 
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U(N) x U(N)/U(N) 


WZNW term 
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CII (chiral symplectic) 
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U(N)/0(N) 


Z 2 term 
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CI (BdG) 


+1 


-1 
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Sp(2N) x Sp(2N)/Sp(2N) 


WZNW term 
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C (BdG) 





-1 





Sp(2N)/U(N) 


9 term 





Dili (BdG) 


-1 


+1 
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O(N) x 0(N)/0(N) 


WZNW term 
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D (BdG) 





+1 





0(2N)/U(N) 


9 term 






be achieved at the lattice level if the fermion doubling 
problem has been overcome, as is the case with the sur- 
face states of three-dimensional Z-topological band insu- 
lators. 

The Z 2 -topological term in symmetry class All was 
derived in the context of disotdered graphene with long- 
range correlated disorderEaEil or two-dimensional sur- 
faces— of three-dimensional Z 2 topological band insula- 
tors .y 

LeClair and Bernard have extended the RMT classifi- 
cation by demanding that all perturbations to the two- 
dimensional Dirac Hamiltonian with N f flavors preserve 
the Dirac structured In this way, the ten-fold classifica- 
tion can be refined by discriminating the parity of N f for 
the 3 symmetry classes AIII, Dili, and CI. These 3 sub- 
classes correspond to the fact that the replicated princi- 
pal chiral models (PCMs) whose target space correspond 
to symmetry classes AIII, Dili, and CI, respectively, can 
be augmented by WZNW terms. The realization of any 
of these additional 3 subclasses in a lattice model requires 
overcoming the fermion doubling problem. 

The parity of the flavor number N f of random Dirac 
fcrmions also matters for symmetry classes All and CII. 
The fermionic replicated NLcrMs derived from the ran- 
dom Dirac Hamiltonians in symmetry classes All and 
CII can acquire a Z 2 topological term on account of the 
dimensionality of the Dirac matrices (twice the number 
N f of flavors) that represents the random Dirac Hamil- 
tonian. Deriving these Z 2 topological terms from lattice 
models is not automatic, for the fermion doubling prob- 
lem must be surmounted. 



In this paper, by identifying a disordered fermionic 
model that gives rise to the Z 2 -topological term in sym- 
metry class CII, we complete the derivation for non- 
interacting fermions subject to a weak white-noise cor- 
related random potential of topological or WZNW terms 
in all 10 symmetry classes relevant to two-dimensional 
Anderson localization. The microscopic fermionic model 
is realized by the surface states of a three-dimensional 
Z 2 topological band insulator in symmetry class CII of 
Rcf. [ll]. (Sec Rcf. ^3| for a particular lattice model of a 
three-dimensional Z2 topological insulator in symmetry 
class CII.) 



C. Global phase diagram 

In this paper, we start from the kinetic Hamiltonian 
JC for N f = 2 flavors of Dirac fermions that make up a 
(reducible) 4-dimensional representation of the homoge- 
neous Lorentz group SO(l,2). We then subject JC to a 
static and chiral-symmetric random potential V, i.e., the 
random Dirac Hamiltonian % = JC + V must anticom- 
mute with a unitary matrix C, {"H, C} = 0, which squares 
to the identity. By imposing the condition that % is in- 
variant under a representation T = of time reversal 
for spinless single-particle states, % belongs to symme- 
try class BDI in the 10- fold classification (see Table |). 
This corresponds to an antiunitary time reversal operator 
whose square equals plus the identity. 

It is also known that such a Hamiltonian % describes 
graphene (see Fig. 0) or the two-dimensional 7r-flux 



4 



E(k) 



(a) 



(b) 











-71 


12 


+71 


12 




E(k) 











-71 


/2 


+71 


/2 




m(x) 



+m 



-m 



(e) 





FIG. 2. (a) Ring along which a spinless electron hops be- 
tween nearest-neighbor sites shown as circles with the uni- 
form real-valued amplitude t. The lattice sites are colored in 
black on one sublattice and white on the other sublattice. (b) 
Electronic dispersion corresponding to (a) after folding the 
Brillouin zone, (c) Ring along which a spinless electron hops 
between nearest-neighbor sites with the dimerized real-valued 
amplitude t ± St. There are two defective sites belonging to 
opposite sublattices at which two strong bonds t + St meet, 
(d) The breaking of translation invariance in (c) has opened a 
gap at the reduced zone boundaries and localized two bound 
states around the two defective sites, (e) A generalization 
of (c) and (d) in three dimensions can be achieved with the 
help of a suitable dimerization of the 7r-fiux phase depicted in 
Fig. ^(b) for spinless electrons. The continuum approxima- 
tion yields a massive 8x8 Dirac equation. Two-dimensional 
defective surfaces normal to the direction x, say, occur when 
the mass changes sign. One mid-gap state is bound to each 
of the two-dimensional defective surfaces. Each midgap state 
obeys a 4 x 4 two-dimensional massless Dirac equation as de- 
picted by a Dirac cone. The two mid-gap states have opposite 
chiralities. 



phase, in the presence of real-valued, nearest-neighbor, 
spin- independent, random hopping amplitudes when the 
Fermi energy is at the band center and once the long- 
wave-length limit haSr-been taken with respect ta-the dis- 
crete Fermi points. Eiro For the case of graphencj£j static 
random real- valued nearest-neighbor hopping amplitudes 
are induced by neglectingEil the dynamics of phonons rel- 



FIG. 3. The four independent dimerization patterns for the 
real- valued nearest-neighbor hopping amplitudes of a spinless 
electron on the honeycomb lattice that preserve the sublat- 
tice symmetry and the time-reversal symmetry for a spinless 
particle. The two triangular sublattices of the honeycomb lat- 
tice are distinguished by the coloring of their sites (white or 
black colored circles). Strong and weak bonds are depicted 
by thick and thin lines, respectively. The two independent 
Kekule dimerization patterns (a) and (b) are responsible for 
the opening of a complex-valued gap m in the continuum ap- 
proximation by a 4 x 4 Dirac equation. The two independent 
columnar dimerization patterns (c) and (d) are responsible for 
the emergence of an axial vector gauge field or, equivalently, 
the complex-valued axial gauge field a in the continuum ap- 
proximation by a 4 x 4 Dirac equation. 



ative to that of the electrons to which they couple. We 
emphasize that it is imperative to treat all channels (see 
Fig. [}]) of disorder compatible with the chiral and time- 
reversal symmetries. 

The first result of this paper is that analytical con- 
tinuation of the real- valued random hopping amplitudes 
to imaginary ones in the aforementioned bipartite lattice 
models yields a random Dirac Hamiltonian that belongs 
to symmetry class CII, as it now turns out to obey the 
time-reversal symmetry (TRS) generated by an operator 
T' = — T /T acting on an isospin-1/2 single-particle state. 
This corresponds to an antiunitary time reversal operator 
whose square equals minus the identity. 

Second, we argue that, this random Dirac Hamiltonian 
captures the (nearly) critical localization properties of 
the surface states of a lattice model that, in the clean 
limit, realizes a three-dimensional Z 2 -topological band 
insulator in symmetry class CII. 

More specifically, we show that the phase diagram de- 
picted in Fig. [| encodes the localization properties of the 
random Dirac Hamiltonian % = K, + V when the chiral- 
symmetric random potential V is assigned the three 
possible independent disorder strengths g Re m , g lm m , g a , 
which arc not irrelevant under the RG. Here we discuss 
a "global phase diagram", depicted in Fig. ^(a), in the 
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space of these three couplings which is projected onto 
the SRem-ftmm P lanc ( with 9 a > = 0). In this phase 
diagram, disordered Dirac fermion systems in all three 
chiral symmetry classes, AIII, CII, and BDI occur in 4 
quadrants, sharing one corner which represents the clean 
Dirac fermion limit. Also realized in the phase diagram 
are the symmetry classes All and D at the boundaries 
separating the three chiral symmetry classes, whereby 
the parametrization of class D turns out to follow from 
analytic continuation of the relevant disorder strength 
that parametrizes class All in the phase diagram. 

The random Dirac Hamiltonian H whose potential V is 
restricted to symmetry class All captures the transport 
properties at long wave lengths of the surface states of a 
disordered three-dimensional Z 2 -topological-band insula- 
tor in symmetry class All (say, Bi-L^Sb^.EJ 

The random Dirac Hamiltonian H whose potential V 
is restricted to symmetry class D captures the trans- 
port properties of the fermionic quasiparticles of a dis- 
ordered two-dimensional chiral p-wave superconductor 
(say, Sr 2 Ru0 4 ) or their counterparts in the random bond 
Ising model at long wave lengths. 

Located in the center of the phase diagram of Fig. ^(a) 
is a vertical dashedJine. There exists a sector of the the- 
ory that decouplesE3 from the random U(l) gauge po- 
tential. This sector is critical along the dashed line in 
Fig. ||(a). We will call the dashed line in Fig. ||(a) a line 
of nearly-critical points to account for the non-critical 
sector that is not depicted in Fig. ^(a). 

It is argued in Sec. |v| that along the dashed line in 
region CII of Fig. ^(a), the transport properties of % are 
also encoded by those of a NLerM on the target manifold 
appropriate for this symmetry class. (Such a possibil- 
ity was also discussed, independently and from a differ- 
ent perspective, in Refs. |64|, and ^5|.) Remark- 
ably, the standard kinetic energy of the NLerM must 
be augmented by a Z 2 -topological term (see Appendix 
|a|). Here, the necessary requirement for the presence of 
the Z 2 -topological term is that the number N f of flavors 
be two times an odd integer n, i.e. N f = In. How- 
ever, any purely two-dimensional non-interacting local 
tight-binding Hamiltonian with Fermi points at the band 
center that breaks the spin-rotation symmetry but pre- 
serves the time-reversal and sublattice symmetries yields 
a Dirac equation with N f = 2n where n is an even integer 
because of the fermion doubling problem. The fermion 
doubling problem for fermions in two dimensions can be 
circumvented by working with fermions localized at the 
two-dimensional boundary of a three-dimensional crystal, 
i.e., with the boundary states of a topological band in- 
sulator in symmetry class CII. It is the nearly-critical lo- 
calization properties of these surface states that are cap- 
tured by the dashed line in region CII of Fig. | Thus, 
we can view the Z 2 -topological term in the NLerM for 
symmetry class CII as the signature of the physics of 
(dc)localization, that arises from the existence of bound- 
ary states in the clean limit, the defining property of 
three-dimensional Z 2 -topological band insulators in sym- 



metry class CII. 

Third, we argue that the initial flow away from the 
apparently unstable nearly-critical line in region CII de- 
picted in Fig. ^(a) is not a crossover flow to the diffusive 
metallic fixed point of the NLerM in symmetry class All 
augmented by a Z 2 topological term. Rather, it is the 
flow depicted in Fig. |(c) that bends back towards the 
nearly-critical plane defined by the dashed line and the 
out-of-plane axis for the coupling g a , as a result of the 
RG flow of the coupling g a , to strong-coupling. This flow 
on sufficiently large length scales along trajectories in 
the three-dimensional coupling space is depicted through 
the two-dimensional cuts presented in Figs, ^(b), ^(c), 
and |(d). The full RG flow along the boundary All, a 
separatrix of the RG flow, was computed numerically in 
Refs. |6^ and ^ owing to the presence of a Z 2 -topological 
term on the target ma nifol d of the NLerM appropriate 
for symmetry class AII.ElEil 

Finally, in the quadrant labeled by BDI, the dashed, 
line also represents a line of nearly-critical points. Oi£3'l£| 
This line of nearly-critical points is stable, without the 
reentrant behavior of the kind mentioned in the preceding 
paragraph. The one-loop RG flow along the boundary 
D, again a separatrix of the RG flow, was computed in 
Refs. |§ 0, and fn|. 

The fact that the quadrant in symmetry class BDI can 
be analytically continued to the quadrant in symmetry 
class CII suggests that one can compute properties of the 
latter phase from the former one. In particular, sets of 
non-perturbative and exact results have been obtained 
for e.g., boundary multifractal exponents for the point 
contaci-cpnductancc on the critical line in symmetry class 
BDI.E30 These results will also apply to the critical line 
in symmetry class CII upon suitable analytical continu- 
ation. 



D. 



Outline 



The rest of the paper is organized as follows: The 
non-interacting random Dirac fermion model is defined 
in Sec. The main result of this section is captured 
by Fig. J We argue in Sees. [II and IV that the gener- 
ating function for the moments of the retarded Green's 
functions for microscopic parameters corresponding to 
the quadrant CII in Fig. [| realizes a replicated fermionic 
or, alternatively, a supersymmetric (SUSY) NLerM aug- 
mented by a Z 2 -topological term. We conclude in Sec. M. 



II. DEFINITIONS AND PHASE DIAGRAM 



We begin in Sec. [I A by defining a non-interacting 
random Dirac Hamiltonian and proceed with a symmetry 
analysis. To identify the axis of the phase diagram in 
Fig. ^ a generating function for the disorder average over 
products of N retarded single-particle Green's functions 
is needed. This is done using the supersymmetric (SUSY) 
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FIG. 4. Global phase diagram for random Dirac fermion defined by Eqs. fl2.12| ), (gJSj), ( g-lj ), ( |2.23| ), and ( p.29| ). (a) Flows of 
the coupling constants close to the clean Dirac point (the origin denoted by an open circle). Along the boundaries D and All, 



the coupling constant g a 
classes D and All require 



ot generated under the RG, so g a , — can be imposed in a consistent way. In fact, symmetry 
g a , — 0. Away from these boundaries, g a i grows under the RG and we have projected the flows 
onto the g a , = plane in the regime where g a , is still small. In the r egion denoted BDI of the phase diagram, there exists a line 
of (nearly) critical points denoted by a dashed line as a result of Eq. ( ^.59b ). This line of (nearly) critical points is perturbatively 
stable under the RG flow (2.63). In the regi on den oted CII of the phase diagram, there exists a line of (nearly) critical points 
denoted by a dashe d line as a result of Eq. (2.621). This line of (nearly) critical poi nts a ppears to be perturbatively unstable 
under the RG flow ( |2.65| ) for small values of g a ,. (b) Infrared flows dictated by Eq. (2.63) close to the clean Dirac point when 
g a , > 0. The slopes of the flows o n the BDI boundaries (?i mm <?R, cm = have changed as compared to the case when g a , — 0. 
(c) Infrared flows dictated by Eq. (2.65) close to the clean Dirac point when g a , > g' + with g'± := g imm i ±p Rcm < and 

9+ > IS- 1- 

The slopes of the flows on the CII boundaries g[ mm gj{, cm = have changed as compared to the case when g a , — 0. Moreover, 
because of the condition g a , > g' + , the RG flows in the quadrant CII are towards the surface defined by the dashed line of 
(nearly) critical points (the g' + axis) and the out-of-plane g a , axis. The plane <? Imm / - g a i with = and g a , > and the 

plan e q p , - g a , with <? Imm / = and g a , > are always unstable under the one- loop flow ( ^.65[ ). (d) Infrared RG flows of 
Eq. (2.65) in the surface defined by the g'_ axis as horizontal axis and the g a , axis as vertical axis of the quadrant CII. 



formalism in Sees. II B and II C . The flows in Fig. |]to or 
away fr om t he nearly-critical line follow once it is shown 
in Sec. IIP that the SUSY generating function defines a 
gl(2N|2N) fc=1 SUSY Thirring model studied in Refs. || 
and (5^. 

A. Definitions 

Common to all the aforementioned microscopic exam- 
ples is the existence of 4 Fermi points at the relevant 
Fermi energy around which linearization in momentum 
space yields the continuum Dirac kinetic energy 




(2.1) 



VxPx + VyPy 


up to a unitary transformation. Here, the momentum 
p = (p x ,p y ) = ip\,P2) is measured relative to the 
Fermi points at the band center. The complex notation 
p = p x — ip y and p — p x + ip y is occasionally used for 
conciseness. The unit 2x2 matrix cr n and the three Pauli 



matrices (a 11 a 2 , <r 3 ) are reserved for the spinor indices of 
SO(l, 2). The unit 2x2 matrix p and the three Pauli 
matrices (p l5 p 2l p 3 ) are reserved for the two-dimensional 
flavor subspace. 

This kinetic energy has two interesting properties. 
First, it anticommutes with the 4x4 unitary and Her- 
mitian matrices 



(2.2) 



Second, the operations on K consisting in the momentum 
inversion p — > p, complex conjugation, and matrix mul- 
tiplication from the left and from the right by the 4x4 
unitary and Hcrmitian matrices 





■= Pa « 


5cr , 


C 1 C 1 


= c x c, = 


+C.C* 


= 1, 


c 2 


:= p 2 d 


5cr n , 


c 2 c\ 


— — C^C^ — 




= 1, 


c. 


■= Po<5 




^3^3 


= C3C3 = 




= 1, 




:= p 1 G 


5^3. 


C4C4 


= C4C4 = 




= 1. 





■= p 3 $ 


5o-i, 




= T 1 T 1 = 


+71 V 


= 1, 


r 2 


:= p i 




r 2 rl 


= T 2 T 2 = 


-r 2 r; 


= 1, 


T 3 


■= Pi 6 






= T 3 % = 


-%v 


= u 




:= P2 6 < 


5o- l5 




= %% = 


-%t: 


= 1, 



all yield K again. For any i,j = 1, • • • 
C,/C(p)C, =-Kip), 



(2.3) 



4, the property 
(2-4) 
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that we will call (abusively) chiral symmetry (chS), is 
compatible with the property 



7} fc*(-p)75 =fC(p), 
that we will call TRS, if and only if 



0. 



(2.5) 



(2.6) 



In this paper, we shall assume that the lattice model 
from which JC(p) emerges imposes the chiral symmetry 
generated by 



C = C X . 

This chiral symmetry commutes with 
T = T X 

and with 

r = %. 



(2.7) 
(2.8) 
(2.9) 



(Observe that T and T' anticommutc. They are not 
compatible.) This leads to two possible forms of TRS, ei- 
ther the one appropriate for particles with integer isospin 
when 



T 1 



-T 



(2.10) 



is imposed as a symmetry, or the one for particles with 
half-integer isospin when 

T' T = -T' (2.11) 

is imposed as a symmetry. Again, the choice between T 
and T' is dictated by the underlying lattice model. 

The most general static random potential that anti- 
commutcs with C is of the form 



v- v 

V = a 1 A 1 + a 2 A 2 + a 3 M 3 + a Q M , 
where the complex-valued 

A l = a 1 — ia'i, 

^2 = a 2 ~ * a 2' 

M 3 = -m 3 - im' 3 , 



(2.12a) 



(2.12b) 



'o- 



represent sources of (static) randomness, i.e., complex- 
valued functions of the space coordinates r € R 2 . (The 
unusual sign conventions is chosen to make contact with 
the notation of Ref. ^.) It yields the random Dirac 
Hamiltonian 

H{r) := (K + V)(r) 

D(r) 
Dt(r) 



ip 1 <£) <r 1 d l — ip 1 (£> a 2 d 2 
p 1 ®a 1 a 1 (r) + p 1 ®a 2 a 2 (r) 

Pi ® o-i a'i(r) +p 2 ®cj 2 a 2( r ) 
Pl ® cr 3 m 3 (r) + Pl ® cr m' (r) 
p 2 ® cr 3 m 3 (r) + p 2 ® CT o TO o( r )- 



(2.12c) 



By construction, Hamiltonian ( [2.12c ) is a member of the 
AIII symmetry class (chiral-unitary symmetry class) of 
Anderson localization in two dimensions. 



When the disorder (2.12b) is restricted to 



A = —ia' € iR, M 3 = -m 3 e R, M = -im e il 



the random Hamiltonian (2.12c) reduces to 



U{r) 



ipi ® o- x d 1 — ipi ® C7 2 9 2 

p 2 ® o-j oi(r) + p 2 ® cr 2 a 2 (r) 

p 1 ® cr 3 m 3 (r) + p 2 ® fi m (r) 



(2.13a) 



(2.13b) 



and hence is invariant under the time-reversal 

TH*(r)T = H(r), T:=p z 2t<T U (2.13c) 



for any realization of the disorder (2.13a). Accordingly, 
this Hamiltonian is a member of the BDI symmetry class 
(chiral-orthogonal symmetry class) in Anderson localiza- 
tion. 

On the other hand, when the disorder ( 2.12b| ) is re- 
stricted to 



A = —ia' € iR, M 3 = -im' 3 e iR, M = m' e 



the random Hamiltonian (2.12c) reduces to 



(2.14a) 



H(r) 



(2.14b) 



ipi <3 cr 1 d 1 — ipi <E> cr 2 d 2 
+ p 2 ®cr 1 oi(r) + p 2 ® a 2 a' 2 (r) 
+ p 2 ®cr 3 m' 3 (r) + p 1 ®<7 m' (r) 

and hence is invariant under the time-reversal 

T'H*{r)T' =H(r), T' := p g> <r 2 , (2.14c) 



for any realization of the disorder ( 2.14a ). Accordingly, 
this Hamiltonian is a member of the CII symmetry class 
(chiral-symplcctic symmetry class) in Anderson localiza- 
tion. 

The BDI case ( |2.13| ) can be derived as the con- 
tinuum limit of a real-valued, nearest-neighbor, spin- 
independent, and random hopping model on a bipartite 
lattice; the honeycomb lattice oLgraphene or the square 
lattice with a 7r-flux phase say.Lj The four-dimensional 
subspace associated to the p's and cr's originates from the 
2-sublattices structure and the 2 non-equivalent Fermi 
points at the band center. The electronic spin plays here 
no role besides an overall degeneracy factor as spin-orbit 
coupling is neglected. In the context of graphenc, the 
random fields a[ and a' 2 are called ripples [sec Fig. ||(c- 
d)],Lj while the random masses m 3 and m are smooth 
bond fluctuations about the Kekule dimerization pattern 
of tie nearest-neighbor hopping amplitude [see Fig. ||(a- 
b)].L3 In the context of the 7r-flux phase, the random 
fields a[ and a' 2 are smooth fluctuations of the nearest- 
neighbor hopping amplitudes about the two wave vectors 
for the two independent staggered dimerization patterns, 
while the random masses m 3 and m are smooth bond 
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TABLE II. Symmet ry co nditions on the static random fields 
in the Hamiltonian ( |2.12| ), For the symmetry classes D and 
All, M 3 M = must hold. 





AIII 


BDI 


CII 


D 


All 


Ax 


a 1 — ia x 


• 1 

—lCb\ 


—zai 








A 2 


a 2 - ia' 2 




~ia' 2 








M 3 


—m 3 — im' 3 


-m 3 


-im' 3 


-m 3 


-im' 3 


M 


m'o — irriQ 


—im,Q 


m' 


—im,Q 


m' 



fluctuations about the two independent columnar dimcr- 
ization pattern.El The CII case (2.14) can be derived as 
the restriction to a two-dimensional boundary of a dis- 
ordered, three-dimensional Z 2 -topological band insulatp; 
in the chiral-symplectic class of Anderson localization.tll 
When the disorder is restricted to 



\ = 0, 



Mo = — ?Tln G 



M = 0, (2.15a) 



we observe that Hamiltonian ( 2.12c ) reduces to 

H = p 1 ®D, (2.15b) 

with 

D = D\ <r x D*ai = -D, (2.15c) 

and can be thought of as a random Hamiltonian belong- 
ing to the symmetry class D (BdG Hamiltonians with 
both time-reversal symmetry and spin-1/2 rotation sym- 
metry broken) in Anderson localization, for H is then 
unitarily equivalent to 



D 
—D 



(2.15d) 



with the unitary transformation (p + ip 2 ) g) a /y/2. 
Finally, when the disorder is restricted to 



A M = 0, 



M, = 0, 



M n 



m G 



(2.16a) 



we observe that Hamiltonian (2.12c) reduces to 

U = Pl ®D (2.16b) 

with 



D 



<T 2 D*a 2 



D, 



(2.16c) 



and can be thought of as a random Hamiltonian belong- 
ing to symmetry class All (a spin-1/2 electron with time- 
reversal symmetry but without spin-rotation symmetry) 
in Anderson localization, f or % c an then be brought to 
the block diagonal form (2.15d ) b y the same unitary 
transformation used to reach (2.15d). 



All four symmetry conditions are summarized in Ta- 
ble [n| The defining conditions on classes D a nd Al l 
can be made slightly more general than in Eqs. ( |2.15a| ) 
and (|2.16a|) as will become clear at the end of Sec. flC 



B. Path integral representation of the 
single-particle Green's function 

In Anderson localization, physical quantities are ex- 
pressed by (products of) the retarded (+if], r\ > 0) and 
advanced (—iff) Green's functions 



g R / A (E) := (E±ir]-ny 



(2.17) 



At the band center E = 0, the retarded and advanced 
Green's functions are related by the chiral symmetry 
through 



c g r (e = o)c = -g a (e = 0). 



(2.18) 



Hence, any arbitrary product of retarded or advanced 
Green's function at the band center equates, up to a sign, 
a product of retarded Green's functions at the band cen- 
ter. From now on we will omit the energy argument of 
the Green's function, bearing in mind that it is always 
fixed to the band center E = 0. 



Because of Eq. ( 2.18| ), it suffices to introduce functional 
integrals for the retarded Green's function defined with 
the help of the SUSY partition function 



Z : — i?p x Z-q , 

Z F := J V[x, x) exp \ i J d 2 r X {irj - U) x ) 
Z B ■= I exp (i / & 2 rUvi-U)i 



(2.19a) 



Here, (x, x) is a P a i r of two independent four-component 
fcrmionic fields, and (£, £) is a pair of four-component 
bosonic fields related by complex conjugation. For any 
?7>0, 



1 



(2.19b) 



holds. The matrix elements of the retarded Green's func- 
tion can be represented as 



iQ R (r,r') = {x(r)x(r')) = {t(r)Z(r')) 



(2.20) 



with (■ • • ) denoting the expectation value taken with the 
partition function Z . 

We now perform the change of integration variables 
from x, x to ■)/> a ^, ?/> a t, ip a in the fermionic sector and 
from £, £ to (3 a \ /3 a , [3^, /3 a in the bosonic sector where 
a = 1, 2 and. 



X ='■ 


1 

v/2tt 




-#2 " 


-#2) > 


X ='■ 


1 


(+iV 2t ^ 






£=: 


1 

v/2tt 


^ 


-02 ~ 




£=: 


1 

v/2^ 


H/? 2t - 


0V (3, 





(2.21) 
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Any corr elatio n function such as the retarded Green's 
function ( ^.20| ) is, under this or any similar change of in- 
tegration variable, to be computed with the SUSY par- 
tition function 

Z=JV[^ t p t 0\ — — ^ — 

x cxp \i J d 2 r x{ip, i>) (it) ~ K) X$, Vol ( 2 -22) 
x exp (i f d 2 r fC9, /?) (i v - H) £(0, /?) 



The message conveyed by Eq. ( 2.22 ) is th at we are free to 
relabel all integration variables in Eq. (2.19a) indepen- 
dently from each other, provided the correct book keep- 
ing with the in tegration variables in the convergent path 
integral (2.19a) is kept. In this conte xt the symbols ~ and 
f on the right-hand side of Eq. ( 2.21 ) are only distinctive 
labels, i.e., here they are not to be confused with c omp lex 
conjugation. The change of integration variable ( 2.21 ) is 
made to bring the effective action to a form identical, to 
that found in Ref. |56| in which important symmetries^ of 
the partition function in the limit r\ = become manifest. 
We also introduce 



a = a 1 — ia 2 = KeA 1 — iKeA 2 , 
a' = a-L — ia 2 = —ImA 1 + i ImA 2 , 
m = m — im 3 = — Im M + i Re M 3 , 
m! = m' — im' 3 = Re M + i Im M 3 , 



(2.23) 



and their complex conjugate a, a', m, and to', in terms 
of which symmetry class BDI is defined by the conditions 

a = 0, tt'eC, me C, to' = 0, (2.24) 
while symmetry class CII is defined by the conditions 

a = 0, a'eC, to = 0, to' G C. (2.25) 

The boundary 

a = a' = 0, Rcto = 0, m' = (2.26) 

between the symmetry classes BDI and AIII belongs to 
symmetry class D. The boundary 



0, 



0, 



Iitito' = 



(2.27) 



between the symmetry classes CII and AIII belongs to 
the symmetry class AIL All four symmetry conditions are 
summarized in Table III. The defining conditions on the 



symmetry classes D and All can be made slightly more 
general than in Eqs. fl2.26|) and (2.27) as will become 
clear at the end of Sec. |I C| . 

With these changes of variables, the partition function 
can be written as 



Z = Z F x Z B at E 



Z F = J V[r\r\i>^] exp(- J & 2 r(c F + C 
Z B =J V[^,(3^J a ,f3 a ] exp 



d 2 r ( £ F 



TABLE III. Symmetry con dition s on the static random fields 
in the generating function (2.28). 





AIII 


BDI 


CII 


D 


All 


a 
















i 




i ■ 7 

a l ~ la 2 


/ ■ / 

a l ~ la 2 


a 1 — ia.2 








in 


m — im^ 


m — im 3 





m m 3 — 





m' 


m'o — im'3 





m'o — im'-?, 





rn'om'^ = 



with the effective action for the fcrmionic part given by 

±-Y.{ rH2d-i(-lTa + a>}^ 

+ ^ [2d-z(-l) a a + a'K (2.28b) 

+ [TO+ (~l) a+ W] ^ at ?/> a 
+ [?7l+ (-l) a+ W] V^V'a} 



U = 



and 

4" = ^ (V 5 V f + - Mi - iMi) > (2.28c) 

and the bosonic part of the effective action given by 

a=l 

+ /3 a t [2d-*(-l) a a + a']/3 a ( 2 .28d) 
+ [m + (-l) a+ W] ^ at /3 a 

+ [m + (-l) a+ W] /3 at ,5 a } 



and 



(2.28c) 



where 2d = di — id% and 2d = d\ + idi. The asymmetry 



between fermions and bosons in C F and £3, a conse- 



quence of the asymmetry betw een the ^'s and /3's on 
the right-hand side of Eq. (2.21), is the price to be paid 
in order to make a GL(2|2) supersymmetry of C F + C B 
explicit, as is shown in Refs. and p6p3 

The N-th moment of the retarded single-particle 
Green's function evaluated at the band center is ob taine d 
by allowing the index a to run from 1 to 2N in Eq. fl2.28| ). 



(2.28a) 



C. Phase diagram 

We now assume that the disorder potentials arc white- 
noise correlated following the Gaussian laws with vanish- 
ing mean and nonvanishing variances 

w(r) = 0, w{r)w{r>) = g w 8 {2) (r - r'). (2.29a) 
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TABLE IV. Symmetry conditions on the (positive) variances 
of the static random fields from Table III . For symmetry class 



D, ffRemSlmm 


= 0. 


For symmetry class All, g Rem 


5lmm' — 


AIII 


BDI 


CII 


D 


All 


^?R,c a 














9lm a 














9Re a ! 


fee a' 


9 Re a ' 








9lm a' 


9lm a' 


9lm a' 








9lic m 


9Hc m 





9n.c m 





9lm m 


9lm m 





9lm m 





5Rem' 





flRem' 





9Rem' 


9\m m' 





9lm m 1 





5lmm' 


Here, S^ 2 \r 


-r') 


is the two-dime 


nsional delta function 



(• ■ • ) represents disorder averaging, 

w € W :={Rea, Ima, Re a', Ima', 

Rem, 1mm, Rem', Imro'}, 



(2.29b) 



and the disorder strengths g w are all positive. We shall 
treat symmetry class BDI defined by 



5Rc a — Sim a — 9r.c m' _ 9lm m' — 

and symmetry class CII defined by 

9Kca = 9 Ima = ffRcm = ftm m = 0- 

Their boundaries 

9lm a' 



and 



=g Rca z 

= 9Re m> 

°=9Rea z 
— 9Rc m 



9lm a 
9lm m' 



9lma - 
9lm m 



9Re a' 
= 9Rem 



9Rca' ~ 
~ 9lm m' 



9iu 



(2.30) 
(2.31) 

(2.32) 
(2.33) 



to symmetry class AIII are in symmetry class D and in 
symmetry class All, respectively. All four symmetry con- 
ditions are summarized in Table IV. The defining condi- 
tions on classes D and All can b e ma de slightly more 
general than in Eqs. (2.32) and ( 2.33| ) as will become 
clear shortly. 

The phase di agram f or th e r ando m Dirac fermio ns dc - 
fined by Eqs. ( pT2| ), (§1]|, ( fl| ), (§]||, and (|2~29| ) 

belongs to the 8-dimensional parameter space 



AIII 



{g w eR\0<g w <oo,weW} (2.34) 



with the origin representing the clean limit. Imposing on 
f2 AIII the constraints summarized in Table |y| yields the 
4-dimensional subspaces 



^BDI C ^AIII' "CII 

and the one-dimensional subspaces 



r2 nTT c r^^jjj, (2.35) 

(2.36) 



^AII C ^AIII 



We arc going to analyze the phase diagram and the pro- 
jected RG flows of its couplings through two-dimensional 
cuts in fi AIII which we will depict with Fig. |[ All those 
cuts belong to the 6-dimensional subspace 



SI- 1 := {g w e n Am \ = g Rca = ,g Tma }. 



(2.37) 



The cuts will involve a plane with the variance of the 
gauge potential a' set to either zero in Fig. [|(a) or a 
nonvanishing value in Figs. ||(b) and ||(c). Wc shall also 
represent the effect of the RG flow to strong coupling of 
the variance of a' on the coupling g'_ :— <? Imm / — <?R em / 
in Fig. |(d). 

To this end, we observe that the quadrant 



9Rem > 0' 



9h 



>0, 



(2.38) 



belongs to symmetry class BDI in Fig. [|(a). The quad- 
rant 



9Rem < 0' 



9h 



<0, 



(2.39) 



in Fig. [|(a) belongs to symmetry class CII as we now 
demonstrate. This is expected from the fact that m' 3 
present in the CII model is the imaginary counterpart of 
mo 3 present in the BDI model. 



We begin with the Lagrangian (2.28b) on which we 
perform the transformation 



2t 



at 



i/' 2 



"■02 



(2.40) 



Under this transformation 

2 



a=l 
2 

£^¥a 



a=l 
2 

E(-i) a+ v at ^ a 

a— 1 a— 1 

2 2 

a— 1 a— 1 

2 2 

Ev> t v.-a^E(- 1 ) a+1 ^ at ^ I 



(2.41) 



while all other terms in Lagrangian ( [2.28b|) re main un- 
changed. We conclude that Lagrangian ( 2.28b ) remains 
unchanged by combining transformation ( 2.40| ) with the 
transformation 



Re m 



i Re i 



Im m 



i Imm'. 



(2.42) 



As the same argument carries thro ugh in the bosonic 
sector by combining transformation (2.42) with 



02 -> -02 



(2.43) 



we conclude that a disorder realization in symmetry clas s 
CII is obtained from the analytical continuation (2.42) 
of a disorder realization in symmetry class BDI when 
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rj = [C-p and £ p are not i nvari ant under the trans- 
formations (|2.40D , ( 2.42 ), and ( 2.43|) 1. Upon disorder av- 
eraging, the analytical continuation (2.42) amounts to 
mapping the CII quadrant 



ffRem' > 0' 



Slmm' > 



(2.44) 



one-to-one into the quadrant ( [2.39 ) through the mapping 



5r c 



-9rc 



(Jin 



.91,: 



(2.45) 



that relates the positive variances <7R em / and 3i mm / in 
symmetry class CII to the negative variances <? Rcm and 
Simm- The remaining quadrants in Fig. 0(a) 



and 



< .9r c 



> .9r c 



> .9r 



-.9r c 



-9ir, 



< g lr 



(2.46) 



(2.47) 



belong to symmetry class AIII as their corresponding dis- 
order potential p 2 <8> (co TO o + °3 m 3) an d Pi <8> (co m o — 
cr 3 m 3 ) are not invariant under neither the time-reversal 
operation T nor the time-reversal operation T'. 
The one-dimensional boundary 



= S-Rem) < .91,: 
of the BDI quadrant, 

< 2Rom> < g ln 



(2.48) 



(2.49) 



belongs to symmetry class D accor ding to Eq. fl2.32j ). The 
onc-dimcnsional boundary 



< 5r 



= 9lr 



(2.50) 



of the CII quadrant (2A4) belongs to symmetry class All 
according to Eq. ( 2.33[) . The one-dimensional boundaries 



< 9r c 



= 9lr 



(2.51) 



and 



However, the signs of the random fields Rem, Imm, 
Rem', and Imm' are innocuous after disorder averaging, 
for these fields are Gau ssian distributed with a vanishing 
mean according to Eq. ( 2.29 ). Hence, a mirror symmetry 
along the vertical axis in Fig. ||(a) must hold. 

The RG flows along the boundaries D and All are 
known and shown in Fig. ^(a). In symmetr y c lass D, the 
RG flow is to the clean Dirac limit (see Refs.]6^, [70| [n], |^ 
and [77]), while the RG flow is to the metallic fixed-point 
in symmetry class All (see Refs. 51, ^(| and |67|) .Ota The 
random vector potentials a 1 — ia[ and a 2 — ia' 2 . are not 
generated under the RG on the boundaries D and AIL 

The RG flows away from the boundaries D sho wn in 
Fig. ^(a) are consistent with the fact that the line (2.52) 
is a stable line of nearly-critical points in the BDI quad- 
rant. As we show below, they also follow from a one-loop 
stability analysis summarized in Fig. f|(b). The RG flows 
away from the boundaries All shown in Fig. ^(a) are a 
more subtle matter. They are drawn to be c onsistent 
with the fact that the nearly-critical line (|2.53| ) appears 
to be unstable in the CII quadrant of Fig. ^(a) when the 
approximation g a , m is used. However, as we show 
below, relaxing this approximation and allowing the RG 
flow to reach length scales such that g a , becomes suffi- 
ciently large changes the flow depicted in Fig. ||(a) to 
that depicted in Fig. |](c). This change is a consequence 
of the flow depicted in Fig. |(d). 



D. The plane R B g M = g Kc 



9ln 



and g a , > 



Consider the line ( 2.53 ) in Fig. ^(a). By combinin g the 
results of Refs. p3q and [56| with the results of Sec. II C , 



we are going to show that this line is a line of nearly- 
critical points. To this end, we shall assume that rotation 
symmetry is preserved at the statistical level. This means 
that we can assume 



9 Re a — 9ln 



9ai 



fRo 



9h 



9a' 



(2.55) 







fee 



< g lr 



(2.52) 



also belong to symmetry classes D and All, respectively, 
as follows from the mirror symmetry about the line 



9 5m = #Re 



(2.53) 



To derive this mirror symmetry, one observe s, wh en rj = 
0, the invariance of the generating function ( 2.28|) under 
the combined transformations (a = 1,2) 

Re m — > Im m, Im m — > —Re m, 
Rem' — > Imm', Imm' — > —Rem'. 



(2.54) 



1. The plane g M = g Re 
We begin with the plane 

< .9m = 5Rom = ftmmi 



> and g a , > 



0< 9 a ', 



(2.56) 



in Fig. |] along which the generating function for the av- 
erage retarded Green's function, which is nothing but the 
gl(2|2) fc=1 Thirring model stu died i n Refs. ^ and [3(| In- 
deed, by setting rj = in Eq. ( 2.28| ) and integrating over 
the random potentials, one finds the partition function 



gl(2|2) 1 



5 il(2| 2)l - S 



V[ip\ip,ip\ip] exp (-S 



gl(2|2) 1 



d z d ; 



2ni V 2?r 

o a , = -jjt(-i) A Jg(-i) B , 

O m = -J A B / B A (-1) A 



9m q 
2ir 



M 



(2.57a) 
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The action 



and 



S'n 



(2.57b) 



dg 



M 



9m 



dl 



0. 



(2.59b) 



r 2 , 



r 2 ) is the action in Eq. ( 2.2S 



(z = r-y 

without disorder when r\ = 0. The capitalized index A 
1 , • • • ,4 carries a grade which is either for A = 1 , 2 or 1 
for A = 3, 4. It is the grade of the indices A and B that 
enters expressions such as (— ) A or (— ) AB . The, grade 
(1) thus corresponds to the bosons (fermions).l3 We are 
using the summation convention over repeated indices 
A, B = 1, • • • , 4. We also have defined the supercurrents 



Jl 



Bt 



Bt 



(2.57c) 



where A, B = 1, • • • ,4 and ip A , ip A , i/> , and i/j A ^ now 
denote bosons for A = 1,2 and fermions for A = 3,4. 
(By allowing the graded indices A and B to run from 1 
to 4N, we can compute the N-th moment of the retarded 

single-particle Green's function.) 

Observe that the integr ation measure in Eq. @.57a| ) 
and the free action ( 2.57bj ) are both invariant under the 
local chiral GL(2|2) x GL(2|2) transformation 



Observe that the coupling constant < g M does not 
flow (we emphasize that this is a non-perturbative re- 
sult) while the coupling constant g a , flows to strong cou- 
pling even when it is initially zero. This is what is meant 
with the statement that the plane defined by Eq. ( 2.56| ) 
(and its projection onto a half-line) is nearly-critical: it 
is critical (in spite of the flow of the coupling g a , ) for all 
correlation functions of fi elds that are unaffected by the 
flow of g a ,. The half- line ( 2.56| ) in Fig. ||(a) belongs to the 
2-dimcnsional symmetry class BDI in the ten-fold clas- 
sification of Anderson localization (see Refs. |il2|-|T4| and 
Appendix ph. 



2. The plane g M , = g Rcm , = g Imm , > and g a , > 



We continue with the plane 



5m ' = 5r c 



5i„ 



> 0, g a , > 0, 



(2.60) 



At 



Btr-iA 



and 



At 



1A 



Ra^Pb 



(2.58a) 



in Fig. [|. The half line obtained from the projection to 
g a , = of this plane is also a line of nearly-critical points 
that now belongs to the two-dimensional symmetry class 
CII in the ten-fold classification of Anderson localization 



(2.58b) (see Refs. |12J-|1J)- Indeed, the counterpart to Eq. ( |2.57| ) 



is 



for any anti-holomorphic L(z) and holomorphic R(z) in 
the fundamental representation of G L(2|2) . The trans- 
formation law of the currents under (2.58a) and ( 2.58b| ) 



Z 



1(212), 



S 



gl(2|2) x 



d z d ; 



( i^-O — ^ M ' 



C 7 D p-lB 

D > 



^A R 



2iri V 2tt 

A / 1 \A 7B 



2- 



O 



M 



(2.61) 



J a L£ J c d L~£. (2.58c) O a , = J A A (-1) A J B B (-1) B , 

Cm = — ^a Jb(~ 1) A ' 



Hence, the Thirring model (2.57) is invariant under 
the g lobal diago nal sub group of the global transforma- 
tion ( 2.58a ) and ( 2.58b ) defined by choosing 



R = L 



(2.58d) 



in Eqs. (2.58a) and (2.58b) to be independent of space. It 
can be shown that the rj term responsible for the conver- 
gence of the integrals in the bosonic sector that has been 
neglected so far breaks this symmetry down to the sub- 
supergroup OSp(2|2). In fact, the symmetry-breaking 
pattern GL(2|2) — > OSp(2|2) occurs due to supcrficld 
bilinears acquiring an expectation value with the conse- 
quence-of-a diverging density of states (DOS) at the band 
center £JO 

The (infrared) beta functions for the couplings g a , and 
g M haverbeen computed non-perturbatively in Rcf. |5^. 
They areEil 



as follows from the analytical continuation g M — > —g M , 
of Eq. (2.57) or b y exp licit integration over the random 
potentials in Eq. ( [2. 28 ) with 77 = 0, whereby one must 
account for the extra imaginary number multiplying the 
random mass m' for symmetry class CII relative to the 
random mass m for symmetry class BDI in Eq. ( p. 28 ). 
Accordingly, the counterparts of Eq. ( [2.59 ) are 



Pa 



9w 



1 -5m'/ 71 " 



and 



0a 



dl 



1 



5m 



t V 1 +5 M / 7r 



(2.59a) 



where one must impose the condition 

< g M , < tt 
to avoid the pole in the beta function for g a , . 



(2.62a) 



(2.62b) 



(2.62c) 
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E. Conjectured RG flows in Fig. ^ 

We are now going to justify why we have conjectured 
the RG flows depicted in Fig. |j. More precisely, we make 
the following claims. 



• The boundaries D and All in the plane g a , 
RG separatrices. 



are 



• The plane defined by the dashed line in Fig. ||(a) 
and the out-of-plane g a , axis is a stable nearly- 
critical plane in that all RG trajectories from re- 
gion BDI or CII, except the fine-tuned RG flows 
along the separatrix D and All, reach this plane 
asymptotically in the infrared limit. 

• The rationale that allows us to deduce from one- 
loop flows nonperturbative statements is that the 
anomalous scaling dimension of the operator that 
couples to the 'asymmetry coupling' g_ = g Rcm — 
5imm m the quadrant BDI, or to g'_ = .9R C m' ~ 
5imm/ m the quadrant CII, is known to all orders 
in g^ (By definition, g_ = on the dashed line in 
Fig-l) 

To substantiate these three claims, we treat first the 
BDI case and then the CII case. 

The stability analysis in region BDJUof Fig. ^ is deter- 
mined by the one-loop RG equations^ 



9+ ~ 9- 

_9l 
4tt' 

(.9+ + 9g') 9- 

47T 



(2.63a) 



condition g + > \g_\ 1 the variance g a , flows to strong cou- 
pling and the RG flows follow three-dimensional trajec- 
tories. We depict them by using a two-dimensional pro- 
jection in Fig. |^(b). T he perturbative flows in the region 
BDI from Eq. ( 2.63a ) afterprojection to the ff Rcm -.9i mm 
plane are depicted in Fig. Q(b) for g a , ^= 0. These flows 
depict the instability of the BDI boundaries g Rcm > 0, 

fim m = and .9r c m = °> Sim m ^ to H infrared flow to- 
wards the nearly-critical plane g_ = 0.EJ It can be shown 
by adapting nonperturbative results from Re f. [55|t hat the 
beta function for the coupling g_ in Eq. ( p. 6^ ) holds 
to all orders in g a , and to linear order in g_.l£3 Hence, 
we conjecture that the infrared flows are from the BDI 
boundaries to the nearly-critical plane spanned by the 
dashed line and the out-of-plane axis g a , in Fig. |](b) for 
the entire quadrant BDI. 

The stability analysis of the region CILof Fig. [| is 
determined by the one-loop RG equations^ 



A, 



9'1 



9- 



where 



9± : = Sir 



4tt ' 
4tt' 

(g+ -gq')g- 
47T 



9+ > \9 



(2.65a) 



(2.65b) 



where 



These one-loop flows must respect the conditions g a , > 
and g' + > \g'_ \ in order to represent the effects of disorder 
on the underlying microscopic Dirac Hamiltonian and are 
valid in the close vicinity to the clean Dirac point g a , = 
g' + = g'_ = denoted by an empty circle in Fig. ||(a). In 
the regime g , = and g' + ,g'_ -C 1, the line defined by 
any one of the two boundaries All from Fig. |](a) becomes 
the separatrix of the Kosterlitz-Thouless flow 



9± 



Sir 



±5Re 



> 



5-1 



(2.63b) 



These one-loop flows must respect the conditions g a , > 
and g + > \g_\ and i5 , R, C m'".9imm' or der to represent the 
effects of disorder on the underlying microscopic Dirac 
Hamiltonian and are valid in the close vicinity to the 
clean Dirac point g a , = g + = g_ = denoted by an 
empty circle in Fig. ||(a). In the regime g a , = and 
9+i .9- 1, the line defined by any one of the two 
boundaries D from Fig. ^(a) becomes the separatrix of 
a Kosterlitz-Thouless flow 



Pa 



9t_ 

"4tt' 
9+9- 

4-7T 



(2.64) 



(a), we plotted the Kosterlitz-Th oules s 



In Fi 

flows ( |2.64| ) which accurately capture the flows (2.63) 
when g , w 0. However, in the region BDI defined by the 



"4tt' 

g'+g'- 

47T 



(2.66) 



In Fig^ 0(a), we plotted the Kosterlitz-Th oules s 

flows ( |2.66| ) which accurately capture the flows (2.65) 
when g , m 0. However, in the region CII defined by 
the condition g' + > \g'_\, the variance g a , flows to strong 
coupling and the RG flows follow three-dimensional tra- 
jectories. We depict them by using two-dimensional pro- 
jections in Figs. |(c) and^ (d). T he perturbative flows in 
the region CII from Eq. ( 2.65a ) after projection to the 
SRem'-Simm' P lane are depicted in Fig. |(c) when g a , is 
large. These flows show the instability of the CII bound- 
aries flRem' >Pl 9l mm > = and .9 R cm' = °. 9jmm> ^ 

to any g a , > 0.lj Moreover, these flows also show the in- 
frared flow towards the nearly-critical plane g'_ = due 
to a reversal in the direction along the g'_ axis of the 
infrared flows caused by the growth of g a , as is depicted 
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in Fig. |J(d). It can be shown by adapting nonperturba- 
tive results from Rcf. |5^ that the change in the sign of 
the beta function for the coupling, g'_ holds to all orders 
in g a , and to linear order in <?VEj Hence, we conjecture 
that the infrared flows emerging from the CII boundaries 
continue to the nearly-critical plane g'_ = in Fig. |](c) 
for the entire quadrant CII. 



III. 



PROJECTED THIRRING MODEL 



We now proceed by discussing the dashed line in Fig. |J, 

9Rcm = ftmm an( l JRem' = 9lmm" 

If we are only interested in correlation functions that 
arc not affected by the flow (to strong coupling) of g a , , 
we can set g a , = in Sec. [□]. This is because, along the 
dashed line, the coupling c/ a , turns outfij to never feed 
into the RG equations for the remaing two couplings, 

SRcm and fflmm; 0r 5r m' and 9\ m m ' 

A mathematically consistent way to achieve this is 
to replace the affine Lie supcralgcbra gl(2|2) 1 by its 
affine Lie subsup cralgc bra ps l(2|2 ) 1 i.e., the gl(2|2) 1 
Thirring models (2.57) and (2.61) are combined into the 
psl(2|2) 1 Thirring model defined by 



~ 1(2 | 2)i = / V[^\if},ip\ip)exp (- 
dzdz g M 



S 



gl(2|2) 



S il(2\2) 1 - S ° 



2ni 2?r 



M • 



(3.1a) 



Ou = -J a B ^b A (-1) A , 
subject to the psl(2|2) 1 constraints 

o = j A A R A = j A A (-y 

and 



= J A 



Ja 



(3.1b) 



(3.1c) 



along the now critical line g M G R. The constraint ( 3.1b| ) 
justifies setting g a , = 0. The sign of the variance g M dis- 
tinguishes symmetry class BDI (g M > 0) from symmetry 
class CII (<7 M < 0). The graded index A runs from 1 
to 4N when dealing with the N-th moment of the single- 
particle Green's function. 



IV. 



RELATIONSHIP TO A NLcrM 



So far, we have relied on a description of the global 
phase diagram and, in particular, of the vertical dashed 
line of nearly-critical points in region CII of Fig. | that 
makes explicit the Dirac structure underlying the clean 
limit of the theory. In this section, we seek an alternative 
description of this line, in particular far away from the 
clean Dirac limit. 

To this end, we first observe that we can derive a repli- 
cated NLcrM by integrating out replicated Dirac fcrmions 



in favor of Goldstone modes as is done in Appendix [A|. 
We find a replicated NLcrM augmented by a term of topo- 
logical origin, the 9 term at 9 — ir. The same calculation 
also applies to the supersymmetric formulation of the dis- 
ordered system, yielding a 9 term at 9 = ir for the NLcrM 



defined on the SUSY target manifold given in Eq. (4.4) 
below. 

Without the 9 term at 9 = it, this NLcrM was al- 
ready derived starting from a different microscopic model 
within the chiral symplectic symmetry class CII by Gade 
in Ref. This NLcrM has two coupling constants t M , 
and t a , that are positive numbers, in addition to the topo- 
logical coupling 9 = 7r. The labels of these couplings are 
chosen to convey the fact that t M , does not flow (Ref. |8^) 
whereas t a , does flow away from its value at the Gaus- 
sian fixed point (Ref. ^), by a nalog y to the flow of the 
couplings g M , and g a , in Eq. ( [2.61 ), respectively. The 
topological coupling does not flow, for it can only take 
discrete values. 

The question we want to address in this section is what 
is the relationship between this NLcrM with a 9 term at 



9 = 7r and the Thirring model defined in Eq. (2.61). We 



are going to argue that they are dual in a sense that 
will become more precise as we proceed. To this end, we 
shall rely on the SUSY descript ion u sed to represent the 
Thirring model defined in Eq. ( 2.61 ). 

We begin by establishing the r elevan t pattern of sym- 
metry breaking. The field theory ( 2.61 ) is a GL(2|2) prin- 
cipal chiral model augmented by a WZNW-.term at level 
k = 1 when the couplings g M , = g a , = O.lJ This means 
that the theory at g M , = g a , = is invariant under the 
symmetry supergroup 



GL(2|2) x GL(2|2). 



(4.1) 



The curre nt-current perturbations for any g M , > in 
Eq. ( |2.61 ) lower this symmetry down to the diagonal su- 
pergroup 



GL(2|2). 



(4.2) 



In turn, the symmetry GL(2|2) can be further reduced if 
fermion bilinears acquire an expectation value, as must 
be the case if the global DOS is nonvanishing at the band 
center due to the disorder. This is in fact what happens 
if the analysis of Refs. [sB] and [5€ along the nearly-critical 
line in the BDI quadrant of Fig. | is repeated for the case 
at hand, with the remaining residual symmetry being 



OSp(2|2). 



(4.3) 



The Goldstone modes associated with this pattern of 
symmetry breaking generate the supermanifold 



GL(2|2)/OS P (2|2), 



(4.4) 



which is nothing but the SUSY target space for a NLcrM 
model in symmetry class CII (see Ref. [l^ and Ap- 
pendix ^ of this paper). The critical vertical dashed 
line in quadrant CII of Fig. |] arises from removing the 
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sector GL(1; R) x U(l) from the field theory 
ensuing projected field theory is given by Eq. (|3.1 



corresponding operation on the target space (4.4) 
NLcrM for symmetry class CII yields the manifol 



The 
). The 
the 



= 1/k 



PSL(2|2)/OSp(2|2) 



PSL(2|2)/SU(2|1) 
U(2|2)/[U(l)xU(2|l)] 
CP 2 ' 1 . (4.5) 



We have used here the isomorphism between OSp(2|2) 
and SU(2|1). By setting all fermionic coordinates to zero 
on this SUSY manifold, one obtains the bosonic subman- 
ifold given by 

Boson-Boson (BB) Fermion-Fermion (FF) 
(non-compact) (compact) 
SU*(2)/S P (2) x SU(2)/SO(2). (4.6) 

(The definition of the group U*(2) is given in Ap- 
pendix ||.) We close this symmetry analysis by recall- 
ing£j that the second h omo topy group of the compact 
part of the submanifold (fOt) is not trivial and given by 



7r 2 [SU(2)/SO(2)] =Z 2 . 



(4.7) 



We are now going to argue that, under certain natural 
assumptions detailed below, the vertical dashed line of 
nearly-critical points in region CII of Fig. [| is described 
by a NLcrM with a 9 term at 9 = tt on the CP 2 ' 1 target 
space [Eq. (gj)]. 

To understand what could prevent the identification of 
the vertical dashed line of nearly-critical points as realiz- 
ing the CP 2 ' 1 NLcrM with 9 term at 9 = n, we are first 
going to review the connection between the 0(3) NLcrM 
with the 9 term at 9 = tt and the SU(2) 1 WZNW fold 
theory perturbed by the current-current interaction. l3 

The 0(3) NLcrM with 6»-tcrm at 9 = tt captures the 
low-energy and long-wave-length excitations of antifcr- 
romagnetic spin-1/2 Heisenbcrg spin chains. This field 
theory is related to the SU(2) 1 WZNW field theory by 
perturbing the latter with a symmetry-breaking poten- 
tial (coupling constant h) , which has the effect of chang- 
ing the target manifold of the principal chiral model, 
at h = 0, to that of the NLcrM, at h = oo. (See Fig. 
||.) When the WZNW model is near its weakly-coupled 
ultra-violet (UV) Gaussian fixed point, the flow of the 
coupling h away from this Gaussian fixed point is the 
strongest and brings the theory into the vicinity of the 
weakly coupled (UV, Gaussian) fixed point of the 0(3) 
NLcrM augmented by a 9 term at 9 = tt. In the vicin- 
ity of the SU(2) 1 WZNW critical point, the symmetry- 
breaking potential (coupling h) reduces to the marginally 
irrelevant current-current interaction up to more irrele- 
vant interactions (some discrete symmetries must here be 
invoked). When the coupling constant A of the SU(2) 
principal chiral model augmented by the level k = 1 
WZNW term is close to its critical value A = 1/k = 1, 
the symmetry-breaking potential generates RG flows that 
are close to those of 0(3) NLcrM with a 9 term at 9 = tt. 



(a) 




> PSL(2|2) 



= 1/k 



oo (b) 




FIG. 5. (Color online) (a) Phase diagram for the SU(2) prin- 
cipal chiral model with the coupling constant A SU ( 2 ) that is 
(i) augmented by the WZNW term at level k — 1 and (ii) 
perturbed by a symmetry-breaking potential with coupling 
h . The critical WZNW theory is at A SU ( 2 ) = 1/k and is rep- 
resented on the upper-left corner of the phase diagram by a 
filled circle. The flow along the upper boundary of the phase 
diagram is that of the marginally irrelevant current-current 
perturbation. The lower-left corner of the phase diagram is 
the Gaussian fixed point of the SU(2) principal chiral model 
augmented by a WZNW term at level k = 1, it is depicted by 
a filled hexagon. The lower-right corner of the phase diagram 
is the Gaussian fixed point of the 0(3) NLcrM with 8 term 
at 6 = n, it is depicted by a filled square, (b) Same as in 
panel (a) except for the replacement of SU(2) by PSL(2|2) 
and of SU(2)/U(1) by CP 2 ' 1 under the assumption that there 
is no more relevant perturbations than the exactly marginal 
current-current perturbation at the upper-lpft-«orner of the 
phase diagram, fflhe left vertical boundaryE3E3 and the di- 
agonal boundaryLa are now lines of critical points (colored in 
red). The diagonal boundary that connects the upper- left to 
the lower-right corner is a line of critical points that is argued 
to have a dual representation in terms of a Thirring model on 
the one hand or a NLcrM with a 6 term at 9 = n on the other 
hand. 



When the coupling constant of the SU(2) principal chi- 
ral model augmented by the level k = 1 WZNW term 
is small, the symmetry-breaking potential generates RG 
flows that drive the theory very close to the weakly cou- 
pled (Gaussian) fixed point of the 0(3) NLcrM with a 9 
term at 9 = tt. The envelope of all these RG flows can be 
thought of as the RG flow from the Gaussian fixed point 
of the 0(3) NLcrM with a 9 term at 9 = tt to the SU(2) X 
WZNW critical point. 

The same argument can also be used to relate the prin- 
cipal chiral models defined on the groups SU(7V) and 
SO(M + N), augmented by a k = 1 WZNW term, to 
the NLcrM with the target manifold SU(7V)/S0(.ZV) and 
SO(M + N)/SO(M) x SO(N), respectively, when aug- 
mented by a 9 term at 9 = tt. This argument is confirmed 
by exact results obtained from Bcthc-Ansatz integrabil- 
ity for these NLctMsM^ 

On the other hand, when the level is larger than one 
(for example, k > 1 arises from a fine-tuned half-integer 
spin chain with spin larger than 1/2), the symmetry- 
breaking potential permits (on symmetry grounds) the 
appearance of terms more relevant than the current- 
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^OSp(4«|4n) 1 




(b) 



FIG. 6. (a) Phase diagram for the SU(2) principal chiral 
model with the coupling constant A SU ( 2 ) that is (i) augmented 
by the WZNW term at level k > 1 and (ii) perturbed by a 
symmetry-breaking potential with coupling h . The critical 
WZNW theory is at A SU ( 2 ) = 1/k and is represented in the 
phase diagram by a filled circle. The Gaussian fixed point 
of the SU(2) principal chiral model augmented by a WZNW 
term at level k > 1 is depicted by a filled hexagon. The Gaus- 
sian fixed point of the 0(3) NLcrM with 8 term at 9 — tt is 
depicted by a filled square. The fact that there are operators 
more relevant than the current-current interaction induced by 
the symmetry-breaking potential is indicated by the presence 
of a third axis in coupling space. This third axis quantifies 
the running of the current-current coupling constant g cc that 
is marginally irrelevant. The critical point of the 0(3) NLcrM 
augmented by a 9 term at 9 — n is depicted by an open circle, 
(b) Counterpart to Fig. ^|(b) for the case of the WZNW model 
on OSp(4n|4n) at level k = 1 perturbed by a symmetry- 
breaking potential that projects this WZNW model to the 
NLcrM model in the symmetry class AIL 



for symmetry class All, in complete analogy with the 
case of the projection discussed above from the WZNW 
model to the NLcrM in symmetry class CII. For the case 
of symmetry class All, consider the WZNW model on 
OSp(4n|4n) at level k = 1. In this case, the coupling 
constant of the principal chiral model flows away from 
the WZNW fixed point down towards the weakly coupled 
WZNW model [see Fig. ||(b)]. Now, we project again to 
the NLcrM model in the symmetry class All with the 
help of the corresponding symmetry-breaking potential 
(with coupling constant h). When this is done for the 
weakly coupled WZNW model, this yields the weakly 
coupled NLcrM in class All, the Wess-Zumino term turn- 
ing into a thcta term at 9 = tt on the All target space [see 
Fig. |6|(b)] . On the other hand, the most relevant operator 
in the vicinity of the fixed point of the WZNW model on 
OSp(4n|4ri) at level k = 1 which has the symmetries of 
the symmetry breaking potential is the current-current 
interaction between the Noether currents. This opera- 
tor is marginally relevant. Thus, the RG flow emerging 
from the unstable fixed point of the WZNW model on 
OSp(4n|4n) at level k = 1 ends up in the infrared at 
the weakly coupled NLcrM in symmetry class All [sec 
Fig. B(b)]. This is one way of understanding that the 
NLcrM in class All with the Z 2 term always flows to weak 
coupling (as discussed in Refs. ^ and 67), for it simply 
inherits this feature from the RG flow of the underlying 
WZNW model. 



current interactionsJn the vicinity of the SU(2) fc>1 
WZNW fixed point.E2l Correspondingly, the flow of the 
0(3) NLcrM with a 9 term at 9 = tt will not reach the 
SU(2) fe>1 WZNW critical point, but reaches a different, 
intermediate fixed point describing the critical behav- 
ior of the 0(3) NLcrM with the 9 term at 9 = tt [see 
Fig. |(a)] . This also happens in the OSp(2|2)/GL(l|l) 
NLcrM with theta term at 6 = tt describing tha-critical 
behavior of the spin-quantum-Hall transition.r 9 lr 2 l 

After these preliminary comments, we proceed to the 
case of interest with PSL(2|2) symmetry. If we assume 
that no relevant or marginal interactions other than the 
marginal current-current interaction are allowed at the 
PSL(2|2) WZNW critical point at level k = 1 when the 
PSL(2|2) x PSL(2|2) symmetry of the WZNW model is 
lowered to its diagonal PSL(2|2) symmetry upon intro- 
duction of the symmetry breaking potential (a natural 
assumption for the level k = 1 case under considera- 
tion), then we obtain with Fig. ||(b) the desired relation 
between the PSL(2|2) WZNW theory perturbed by the 
current-current interaction and the CP 2 ' 1 principal chi- 
ral model with 9 term at 9 = ttx3 By analogy with the 
SU(2) fc>1 WZNW critical point, we do not expect this 
assumption to be fulfilled when the level \k\ > 1. 

A similar projection from the WZNW model onto the 
NLcrM with theta term at 9 = tt can also be implemented 



In summary, based on this reasoning we argue that 
the line of nearly-critical points in region CII of Fig. ^ 
(the vertical dashed line in region CII of Fig. |]) has 
two descriptions; one in terms of the PSL(2|2) WZNW 
model perturbed by current-current interactions, and one 
in terms of the CP 2 ' 1 NLcrM at 9 = tt (Z 2 topological 
term). These descriptions are dual to each other in the 
sense that, in the vicinity of the origin of our global phase 
diagram in Fig. |the PSL(2|2) WZNW model is weakly 
perturbed, whereas the CP 2 ' 1 NLcrM is strongly interact- 
ing. On the other hand, for large values of the coupling 
constant g M , of the current-current interaction about the 
Dirac point, a measure of the distance downwards along 
the dotted line away from the clean Dirac point at the 
center of Fig. ^, the resulting Thirring model is strongly 
interacting whereas the CP 2 ' 1 NLcrM is weakly interact- 
ing. We recall that, because the coupling constant g M , 
is exactly marginal, and £fi is the coupling constant of 
the corresponding NLcrM,E3 it is possible to continuously 
interpolate between these two limits by tuning g M , . (The 
possibility of such a duality was also discussed, indepen- 
dently and from a different perspective, in Ref. |6^, |6^, 
and |6|.) 
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DISCUSSION 



Z 2 topological term in the symmetry class CII 
of two-dimensional Anderson localization 



A systematic study of random Dirac fermions in d- 
dimensional space provides a road-map to uncovering 
universal properties of Anderson localization. This is so 
because random Dirac fermions build a bridge between 
models for Anderson localization that are defined on lat- 
tices - and thus are non-universal - and effective field 
theories (NLcrMs) that solely depend on the underlying 
symmetries and dimensionality of space - and as such are 
universal. 

In one-dimensional space, Dirac fermions generically 
emerge after linearization of the energy dispersion around 
the Fermi energy in the clean limit. The effects of 
weak static disorder are then elegantly encoded by a de- 
scription of quasi-one-dimensional quantum transport in 
terms of, diffusive processes on non-compact symmetric 
spaces.Ota This long-wave length description is suffi- 
ciently fine to account for non-perturbative effects such 

ek 
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as parity effects in the numbers of propagating chanr 
in the chiral symmetry classes AIII, CII, and BDI. 
A parity effect can also be derived for the sympl ectic 
symmetry class All in quasi-one dimension.l 101 H 102 l Al- 
though the latter parity effect is not generic in quasi- 
onc-dimcnsional space because of the fermion-doubling 
obstruction, it is generic on one-dimensional boundaries 
of two-dimensional Z 2 -topological band insulators 
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Dirac fermions are the exception rather than the rule 
in band theory when the dimensionality of space d is 
larger than one. Fine-tuning between the lattice and the 
hopping amplitudes is needed to select a linear energy 
dispersion. There is a parallel to this fact in the context 
of Anderson localization. 

For example, in two-dimensional space, the symmetries 
respected by the static disorder do not enforce, on their 
own, the presence of WZNW or Z 2 -topological terms in 
the NLcrM effective long-wave length description of the 
physics of localization. _ 

Ludwig et alS3 (Nersesyan et al.c2) have shown that 
non-perturbative effects can modify the localization 
properties encoded by the two-dimensional NLcrM with 
a WZNW term in symmetry class AIII when studying 
the random Dirac Hamiltonian with N f = 1 (N f > 1) fla- 
vors. Analogous physics can appear in jSymnietry classes 
Dili and CI in two spatial dimensions OEj'EHa However, 
because of the fermion-doubling obstruction, these con- 
ditions cannot be met in purely two-dimensional lattice 
models for Anderson localization. 

On the other hand, they can always be fulfilled on the 
two-dimensional boundaries of three-dimensional topo- 
logical band insulators (itiip-ti are characterized by an in- 
teger topological index) 11X3 

A similar situation holds for the Z 2 -topological terms. 
The number of Dirac flavors N f matters crucially to 
obtain a Z 2 -topological term in symmetry class All as 



shown by Ryu et al. in Ref. |Tj In the present pa- 
per, we have completed the derivation of topological 
terms of two-dimensional NLcrM by constructing the Z 2 - 
topological term for a NLcrM in symmetry class CII as 
a sign ambiguity in the Pfaffian of disordered Majorana 
spinors. Our derivation suggests that this Z 2 -topological 
term cannot arise from two-dimensional local lattice 
models of Anderson localization because of the fermion- 
doubling obstruction, but requires a three-dimensional 
topological band insulator with two-dimensional bound- 
aries. 



B. Global phase diagram at the band center 

The main results of this paper arc summarized in 
Fig. ^. They apply to N f = 2 flavors of random Dirac 
fermions. 

Figure || should be compared with Fig. 9 from Ref. |22| 
that captures the phase diagram for N f = 1 flavors of 
random Dirac fermions or, more precisely, with its pro- 
jection onto the plane g A — in Ref. |2^ (A A = in the 
notation of Ref. ^2|). The phase diagram in Fig. [|is also 
obtained after projecting a three-dimensional flow to a 
two-dimensional subspace of in coupling constant space. 

The three chiral phases AIII, BDI, and CII in Fig. | 
meet at the origin of the phase diagram. This meeting 
point realizes the clean Dirac limit. We showed that an- 
alytical continuation of the disorder at the level of the 
Dirac fermions allows one to move between the BDI and 
CII phases. However, at the microscopic scale of the two- 
dimensional lattice model that realizes the BDI phase, 
this analytical continuation is meaningless. This is yet 
another manifestation of the fermion-doubling obstruc- 
tion. A realization as a local lattice model of the CII 
phase in Fig. || must go through the two-dimensional 
boundary of a three-dimensional topological band insu- 
lator. 

The quadrant BDI in Fig. |] is fairly well understood 
if we assume that the perturbative flows to the nearly- 
critical .dashed line extend all the way to the boundary 
D. BulkEj and boundaryO multifractality and an analytic 
dependence of the conductance on the disorder strength 
g M at the band center E = 0E3 are governed in the ther- 
modynamic limit by their dependence on g M along the 
nearly-critical dashed line. 

The dashed line in region CII of Fig. | is a line of 
nearly-critical points, each of which is captured by a pro- 
jected Thirring model. We have argued that the strong 
coupling regime of this theory is "dual" to a weakly- 
coupled NLcrM augmented by a Z 2 topological term with 
the target space of symmetry class CII. 

We would like to emphasize that the RG flows in the 
CII quadrant of Fig. [|, first away from and then to the 
nearly-critical plane defined by the dashed line and the 
out-of-plane g a , axis, are perturbative in g'± (nonpertur- 
bative in g a i). They are derived under the assumption 
that no relevant or marginal interactions other than the 
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current-current interactions are allowed. The continu- 
ation of these flows to strong coupling is a conjecture. 
At strong coupling, it is tempting to ask if these two- 
parameter flows might be captured by a NLcrM. Evi- 
dently, a NLcrM whose target space is a symmetric tar- 
get space will not do since it would only be character- 
ized by one running coupling constant. A NLcrM on a 
homogeneous but not symmetric target space with two 
independent coupling constants in addition to the Gade 
term would do. (We refer the reader to Rcf. 105 for a 
systematic study of NLcrM on Riemannnian manifolds, 
of which homogeneous and symmetric spaces are special 
examples, as is exp lained in the context of disordered 
systems in Ref. |106| .) We propose that this scenario is 
captured by a NLcrM with the following homogeneous, 
but not symmetric target space (n is an integer, see Ap- 
pendix |j) 



DOS 



DOS 



GL(2n|2ra)/[OSp(n|ra) x OSp(n|n)]. 



(5.1) 



The situ ation here is analogous to the NLcrM discussed in 
Rcf. 106 in the context of the random-bond Ising model in 



two dimensions. The NLcrM on the homogeneous target 
space in Eq. (5.1) has two coupling constants [in addition 



to the coupling constant of the "Gade" term ("projected 
out" in our global phase diagram in Fig. ^), of the kind 
that we previously denoted by g a , in the present paper]. 
The NLcrM on this homogeneous space interpolates be- 
tween the two NLcrMs on the symmetric target spaces 
corresponding to symmetry classes CII and All, which 
are specific limits within the 2-parameter coupling con- 
stant space of the NLcrM on this homogeneous spac e (in a 
manner analogous to the situation discussed in Ref. 106). 
See Appendix E for more details. 



C. Weak breaking of the chiral symmetry in the 
vicinity of the band center 

The effects of a finite Fermi energy E F on the physics 
of localization for the quadrants BDI and CII in Fig. |] 
are dramatic in that, in both cases, a finite E F breaks 
the chiral symmetry chS. 

Turning on a finite Fermi energy E F in the quadrant 
BDI in Fig. || reduces the symmetry—class to AI. All 
states at finite E F arc then localized. El The band cen- 
ter is a quantum critical point separating two insulating 
phases, one defined by E F < and another one denned by 
E F > 0, very much as is the case in the integer quantum 
Hall effect (IQHE) (see Ref. |^ for a review on plateau 
transitions in the IQHE). The global density of states 
v(E F ) diverges as 

»{Ef) ~ exp (-c| In |£ F || 2 / 3 ) (5.2) 

with c a noruajnivcrsal number when E F approaches the 
band center £3 




FIG. 7. (a) Global density of states (DOS) of a two- 
dimensional disorder-free tight-binding model with sublattice 
symmetry and an even number equal to or larger than 2 of 
non-equivalent discrete Fermi points at the band center. The 
dashed line is the global DOS of the corresponding disorder- 
free Dirac fermions in two dimensions, (b) Effect of weak 
disorder for symmetry class All without the Z 2 topological 
term (the surface states of a three-dimensional time-reversal- 
symmetric weak topological insulator), the symmetry class 
relevant to the quadrant CII when perturbed by a finite chem- 
ical potential. Th e b and center is a critical energy at which, 
according to Eq. (5.3), the global DOS diverges. This critical 
energy separates two metallic phases. 



Turning on a finite Fermi energy E F in the quadrant 
CII in Fig. |] reduces the symmetry class to All but with- 
out the Z 2 topological term. Indeed, the random Dirac 
Hamiltonian at a finite chemical potential has now two 
flavors that are coupled by the disorder. This corre- 
sponds to two Dirac cones in any underlying microscopic 
model that arc gcncrically coupled by the disorder. The 
global density of states v(E F ) is again diverging accord- 
ing to the law 



E, 



exp 



(-A 



ln|£ F || 2/3 



(5.3) 



with c' a non-universal number when E F approaches the 
band center.E3 The state at the band center is critical. 
(The robustness to strong disorder of the critical behavior 
of the band center in the chiral classes is well documented 
in quasi-one and two dimensions J123) However, contrary 
to the quadrant BDI in Fig [|, the band center is not 
any more a quantum critical point separating two insu- 
lating phases. Indeed, the localized nature as a function 
of the chemical potential of these two-dimensional states 
is that of the surface states of a three-dimensional time- 
reversal-symmctric weak topological insulator. The issue 
of Anderson localization function of the chemical 
potentia l for suc h su rface states was recently discussed 
in Refs. |l08| and 109. According to the numerical study 
in Ref. |109| (corresponding to the ca se of a mean value 



m = of the random mass m in Rcf. 109 ) , these surface 
states remain extended (metallic) in the presence of dis- 
order even though the characteristic energy at which the 
upturn of the diverging global DOS becomes sizable rela- 
tive to the clean DOS shown in Fig. f?](a) is exponentially 
small for weak disorder £3 
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Appendix A: The sign ambiguity of a Pfaffian 

In this section, we are going to argue that the dashed 
line in the phase diagram of Fig. [| that belongs to the 
chiral symplectic class CII has the particularity that, 
within the fermionic replica NLcrM representation, there 
appears a Z 2 -topological term in addition to the standard 
kinetic energy. 

To this end, it will be useful to enlarge the dimension- 
ality of the representation of the Dirac Hamiltonian by a 
factor of 2 in order to treat the isospin-1/2 TRS. To avoid 
ambiguities, we will use the Greek letters for the Pauli 
matrices acting on the three relevant two-dimensional 
subspaces - p in flavor subspace, a in Lorentz subspace, 
and t in the time-reversal subspace to be introduced be- 
low - as subindices to specify the chosen representations. 
In this section we use indices x, y, z, instead of 1,2, 3, in 
the a and r subspaces. For example, we shall denote the 



Dirac Hamiltonian (2.12) when Eq. (2.14) holds by 



U 



/>.<T 



D c 
Dl 



D a := a x {-id x + A x ) + a y (-id y + A y 
+ a z M z + a M , 



(Ala) 



where A = —ia' € iR, M z = -im' z £ iR, M = mj, € 
R, and with the simultaneous chS 



(p„ ® <r )U p Jp z ® <r ) 



-H 



p,a 



and isospin-1/2 TRS 



(ip ® <J y )n^ a {-ip ®a y ) = H pa . 



(Alb) 



(Ale) 



1. Fermionic functional integral representation of 
the retarded Green's function 

The generating function for the retarded Green's func- 
tion is the partition function 



Z:= J V[x,x] 6x p(-/ d2rC 



Here, we have chosen 



X = (Xi X 2 ) p = (Xit Xu X2f X24.) 



(A2a) 



(A2b) 



to be a 4-componcnt row spinor with Grassmann-valued 
entries. Similarly, 



Xi\ = 

\2, 



Xat 
\X 2 J 



(A2c) 



is a 4-component column spinor with Grassmann-valued 
entries. All 8 Grassmann-valued entries labeled by 
the flavor indices 1 and 2 on which the matrices 
(p , p xl p yl p z ) act and by the Lorentz indices f and J, on 
which the matrices (a , a x , a , a z ) act are independent. 
For the retarded Green's function, 77 > 0. 

It is useful to m ake the TRS (Ale) explicit. To this 



end, following Rcf. [110^ we make the manipulation 

C = - i x(iri - U) Pt<7 x 

= +*x T (*v-u)l <T x T 

= - i X T {-'iPo ® o y )^ - H) Pi<T (-ip ® a y X T ) 



= - 1 m-q - H) p ^ 



(A3a) 



by which we have doubled the number of Grassmann- 
valued entries in $ and ^ through the definitions 



{irj - H) p T a . :=(ir)-H) P!<T 



'0i 



(A3b) 



and 



V2 



1 
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*r *i x l 
( ^ \ 

Xi 

-T 

-xi 

V 4 J r. 



T 

-Xf 



(A3c) 



Here, the subindex r denotes the, by now, explicit time- 
reversal subspace that is spanned by the unit 2x2 ma- 
trix t and the three Pauli matrices (r x , r 2 , r 3 ). Of 
course, the number of independent Grassmann-valued en- 
tries remains unchanged in the representation ( |A3[ ) as the 
TRS ( Ale ) is now represented by the constraint 



tt = tf T (-»p„®T B ®£rA 



(A4) 
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On the other hand, the representation of the chS ( Alb| ) 



is 



t ® a n pTaPz 8 t $ ct = -n PiTilT . 



(A5) 



Instead of Eq. (A4), we seek a representation of the 
TRS in terms of 8-component Grassmann- valued spinors 
obeying the Majorana constraint 



V> = ^ T (-ip ® r (8 cr ). 



(A6) 



This can be achieved by observing that the "square root" 
of T a is given by 



■ T 

-IT T 

1 z-y 1 z-yi 



V2 



(A7) 



Now, we take advantage of the fact that the kernel (A3b) 
commutes with 



With the help of the identity 

i> 2 Dl^ = -^D*^ = ^D T ^ 2 , (A12) 
we arrive at 

C = 2 ^D T ^ 2 - irj (fafa + V^ 2 ) . (A13) 

This presentation of the Lagrangian reveals that, 
upon quantization, ip 1 and ip 2 form a canonical pair of 
fermionic operators. In other words, because of the chS, 
the kinetic part ij>\D T a i\) 2 of the Lagrangian is invariant 
under any global U(2) transformation 

^ -> ^(ctq <g> Ul), ip 2 -> (cr <g> C/ r ) V 2 , (A14) 

where J7 r is a 2 x 2 unitary matrix acting in the r sub- 
space. 



^ ~ Po ® T,_ y <8> <r (A8) 

so that 

£= -»¥(«;- ft), )Tjff tt 

= * T Tj_ y ( lPo ® r ® cr y )(i77 - H) p ^ a T z _ y ^ (A9a) 
= -$(ir)-'H) p ^ a i) 

where 



z-y 



(A9b) 



determines -0 through the Majorana constraint Eq. ( |A6| ) 
that follows because of the isospin-1/2 TRS. In view of 
the Majorana constraint (A6), the isospin-1/2 TRS is 
now equivalent to the global 0(2) invariancc under the 
transformation 

■0 -> (p Q ® CT ® O^), V -> (P ® °0 ® °r ) V> ( A1 °) 

for any 2x2 orthogonal matrix O t acting in the r sub- 
space. 



Finally, it is time to make use of the chS ( Alb ) . By 
making the flavor subspace explicit, 



-irj (V^i + ^ 2 ^ 2 ) > 



(Alia) 



Replicas and disorder averaging 



W o no w assume that a' x , a', 



and 



from 



Eq. (All) are all white- noise distributed with the same 
variance g. In doing so, we limit ourselves to the nearly- 
critical line of region CII in Fig. ||. 
We replicate the Lagrangian N r times, 



2N 

E 

a = l 



[2^ al (~ia p d p + V)^ a2 

ir i (^alV'al + ^2^2) ] ■ 



(A15) 



This Lagrangian is invariant under any global 0(2N r ) 
rotation in the r and replica subspaces. After disorder 
averaging has been performed, we arrive at the interact- 
ing Lagrangian 



2N„ 



£ N =2VdI(-io-,A)d=, 



a=l 

2N r 

f ir]^2 [ d l ifJ v ( d l) T + d Iiv y d a 

3 = 1 

4 - 

,b=l 



(Al6a) 



E (s 



where (a',m' z ,m' £ R) 



D r,a ■= T o 



V) 



V := ia„a[ 



(Allb) 



acts on the two independent 4-component Grassmann- 
valued spinors ip 1 and ip 2 while the spinors ip 1 and ip 2 
obey the Majorana condition (S := r ® <r y ) 



"01 = ^("ffl J, "02 = ^H^J- 



(Allc) 



Here, since tp 1 and 2 are canonically conjugate, we have 
introduced the following notation for any a = 1, • • • , 2N r , 



dl 



S 3 ■= -jd\ad 31 



:= did, 



(A16b) 



It is worth remembering t hat th e replicated "spin" S 3 
in the t- J- like Lagrangian (Al6a) originates from the a 
subspace and not the true electronic spin-1/2. When 
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r/ = and in accordance with the global symmetry ( A14), 
the action is invariant under any global U(2IM r ) rotation 



4 4ui, 



U1U. 



ab ac 



(A17) 



while, in accordance with the global symmetry (AlC), 
any non-zero 77 breaks this symmetry down to the global 
0(2N r ) rotation 



d I -> d b° a b! d a -> °ac d c; °ab°ac = 5 bc, 



(A18) 



where summation over repeated indices is assumed. 



3. Hubbard-Stratonovich transformation 

It is time to introduce auxiliary (Hubbard- 
Stratonovich) fields that decouple the interactions 
among replicas. A possible channel for decoupling is sin- 
glet superconductivity as it is favored by the symmetry 
breaking term 77. Hence for any a, b = 1, • • • ,2N r , we 
introduce the order parameters 



-1 

71 
-1 

1 



0\ b :=-^diia y (4) T 



— 1 Ijt ,)t _ jf jt 



Cab := -^dJi<T y d b 

= —/^ (^af^bj. — ^al^bf) ' 

in terms of which the "exchange term" becomes 

S 3 ■ S b - -n 3 n b = -Ol h O ab , 
and, in turn, the Lagrangian becomes 

2N r 



(A19) 



(A20) 



a = l 



2N, 



(A21) 



a=l 



2N 



E 



,b=l 



The interacting term is then decoupled by the 2N r x 2N r 
Hubbard-Stratonovich field A ab and its complex conju- 
gate A a * b , 



2N. 



C K = 2 J24 H ct M 9 m) d a 
a = l 

2N r 



a=l 



2N„ 



a,b=l 



(A22) 



No approximation has y et been invoked. As the in- 
teracting Lagrangian ( |A22| ) is not readily tractable, we 
shall restrict the path integral to slowly varying bosonic 
degrees of freedom (Nambu-Goldstonc bosons). W e firs t 
look for a diffusive saddle point of the Lagrangian ( |A22j ) . 
In the diffusive regime, the auxiliary field A ab (sponta- 
neously) breaks the global U(2N r ) symmetry, along the 
symmetry breaking "direction" controlled by the symme- 
try breaking term 7/. Thus, the spatially homogeneous 
configuration 



1 0ab 



-*|A |*ab 



(A23) 



should be a representative diffusive saddle point, where 
|A | G K is determined from the self-consistent equation 



In 



A 



2<J 



(A24) 



and A is an ultra-violet cutoff. 

This choice of a saddle point is not exhaustive. Generic 
saddle points can be con struc ted by making use of the 
global U(2N r ) symmetry ( A17 ) of the kinetic energy: 



2N r 

A ab = ]T U ap A 



2N r 



p,q=i 



0pq [/ bq = -7|A |^[/ ap [/ bp (A25) 
P =i 



where U G U(2N r ). However, not all U G U(2N r ) gener- 
ate a new saddle point configuration. If U G 0(2IM r ), A ab 
coincides with the reference c onfig uration A 0ab owing to 
the global 0(2N r ) symmetry (A18). This means that the 
set of saddle points A ab is the coset manifold 



G/H = U(2N r )/0(2N r ), 
whose elements can be parametrized by 



UU 1 



U G U(2N r ). 



(A26a) 



(A26b) 



Note that since UU T is symmetric and unitary, 
U(2N r )/0(2IM r ) is a set of symmetric and unitary ma- 
trices. 

We now include fluctuations around the saddle 



points (A25). Since the longitudinal fluctuations (i.e., 
fluctuations that changes |A |) are gaped, wc shall freeze 
them and only consider the transverse fluctuations 



2N 



A ab (r) = -z |A | Y, U 3p (r)U bp (r) 
P =i 



(A27) 



with U(r) G U(2N r ). With the help of the Nambu rep- 
resentation, the effective Lagrangian becomes 



eft 



2N r 

E 

a,b=l 



%D ab [A] lb , 



where 

7, = (dt <F(-ia y )) a , 



7a = 



d 

ia y (d^ 



(A28a) 



(A28b) 
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are related by the Majorana condition 

7 a = {-i°y®ir yla ) T (A28c) 
with t acting in the Nambu space, and the kernel is 



A£ a (r) 



A ba (r) 

+^ab^^ 



(A28d) 



(We ha ve ab sorbed i\/2rj in a rescaling of A.) Because 



of Eq. (A27) A ab = A ba and thus D ab [A] is Hermi tian 



Observe that the eigenvalues of the kernel ( A28d ) are 
real-valued and the nonvanishing ones come in pairs of 
opposite sign. Indeed, we could have equally well pre- 
sented the effective Lagrangian (A28) as 



Arff = ( d f 



d 1 



K i<y y A(r) \ ( d 
-ia y A\r) -K T ){(dT 

(A29) 

where the kinetic energy K was denned in Eq. (2.1) 
and we use a matrix convention to make explicit the 
Bogoliubov-de-Gcnnes particle-hole symmetry responsi- 
ble for the aforementioned pairing of eigenvalues. 

The effective field theory S cS [A] describing the dynam- 
ics of the slowly varying bosonic field A ab follows from 
integrating out the fermionic fields d) and d in the par- 
tition function, 



-Srff[A] 



V[d\d]cxp(- / d 2 r£ cff 



:(±)v / DetD[A] 
= PfZ>[A]. 



(A30) 



Here, the Pfaffian PfD[A] implements theisospin- 1/2 
TRS through the Majorana condition (A28c) [see also 
Eq. ( Allc )l- A gradient expansion of the exponentiated 
Pfaffian gives the standard kinetic energy of the NLerM 
on the target space U(2N r )/0(2N r ). However, since the 
second homotopy group of G/H = U(2N r )/0(2N r ) is 
non-trivial, 

tt 2 [SU(M)/0(M)] = Z 2 , for M > 2, (A31) 

the NLcrM is allowed to ha ve a t opological term of the Z 2 
type. In other words, Eq. ( A31 ) tells us that the space of 
all field configurations is divided into two sectors that are 
not smoothly connected. Consequently, these two sec- 
tors can be weighted differently in the effective partition 
function. This possibility is enco ded in the ambiguity in 
defining the sign of the Pfaffian ( A30 ), a global property 
of the target manifold G/H = U(2N r )/0(2N r ). In the 
following, we use the same approach as in Ref. to 
show that the ambiguity in defining the sign of the Pfaf- 
fian can be interpreted as the presence of a Z 2 -topological 
term. 



4. Z 2 configurations of the A-field 

In this section, we construct representative A-field con- 
figurations that belong to the two complementary Z 2 - 
topological sectors as defined by the second homotopy 



group (A31). To this end, we introduce the generator 

G o(2N r ) (A32) 



A, := 











2N„ 



of the symme try-b roken group U(2N r ) that leaves the 
saddle-points ( A25 ) invariant. We also define the gener- 
ators X 1 and A 3 through 




(A33) 



2N.-2 



Here, 2N _ 2 is the (N r — 2) x (N r — 2) matrix with in 
all entries. The three matrices A 1: A 2 , and A 3 generate 
an SU(2) algebra. Unlike A 2 , neither \ x nor A 3 leave the 
saddle-points ( A25| ) invariant. Hence, neither X x nor A 3 
belong to the unbroken symmetry group H = 0(2IM r ). 

Let S 2 denote the two-sphere and choose the polar 
— 7r/2 < 9 < +7r/2 and azimuthal < <fi < 2tt angles as 
sphe rical coordinates of S 2 . Following Weinberg et al. in 
we define on S 2 the unitary matrices 



Ref. 



Ill 



U l (0, cj)) := e u ^ /2 e lX ^' 2 e" 11 ^' 2 G U(2N r ) (A34) 

that we label by the integer I G Z. Finally, we define the 
family 



:= -z|A 1 



R, 



= -i|A„ 



in U(2N r )/0(2N r ) where 







'2N 



(A35a) 



cos o + i sm f cos / 
—ism.6sm.ld) 



—ism9smld> 
cos 9 — i sin 9 cos ld> 

(A35b) 

is labeled by the integer I G Z but is independent of the 
replica number N r . 

The Z 2 configurations of the A-field on the two- 
dimensional torus T 2 with the coordinates < x, y < L 
(L is serving as an infrared cutoff) can be obtained from 
the parametrization of the unit sphere S 2 in terms of the 
L-periodic unit vector 



n{x,y) 



r(x,y) 



\r{x,y)\' 

itself given by the L-periodic vector 

/ -sin(27ry/L) 
r(x,y) := — sm(2nx/L) 

\cos(2ttx/L) + cos(2ny / L) 



(A36a) 



(A36b) 



1 



For example, the L-periodic A l=1 (x,y) is obtained from 
Eq. (A35) by replacing R l=1 (9,cj>) with 

R l=1 (x,y) = 



n z + m x 
—in,. 



in. 



(A37) 
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FIG. 8. The energy eigenvalue spectrum in the vicinity 
of th e ban d center for the kernel D[A(t)] from Eqs. ( jA28d ) 
and (A39), is computed numerically as a function of the pa- 
rameter < t < 1 for |A | /A = 1. The field A(t) interpolates 
between Ai when t = and when t = 1. Without loss of 
generality, the replica number N r = 1 was chosen. 



6. Summary 

In summary, after integration over the Majorana 
spinors along the nearly-critical line of region CII in 
Fig. |[ the effective action for the Nambu-Goldstone field 
A, a symmetric and unitary matrix, is given by 



^topolo 
j NLctM 



P[A](-1) 



-S[A] 



(A40a) 



where S [A] is the (fermionic replica version of tlae) action 
for the NLcrM on G/H = U(2N r )/0(2N r ), i.e.,Ell 



S[A]=-L [ d 2 rtr[(9 M At) (^A)] 

l M' J 

+ — f d 2 rtr(At,9 Al A)tr(Aa A1 At) 



(A40b) 



while 



Spectral flow 



We are going to argue numerically that 



B gnPf(D[A,]) = -8gnPf(D[A I+1 ]) J ieZ, (A38) 



n[A] =0,1 (A40c) 

- the Z 2 -topological quantum number of A - reflects 
the ambiguity in defining globally the sign of the Pfaf- 
fian o f Maj orana spinors. Because of the block struc- 
ture ( A35bj ), the topological quantum number n[A] = 
0, 1 is expected to survive the replica limit N r — > 0. 



by looking at the spectral flow of the kernel 



Appendix B: Patterns of symmetry breaking and 
supermanifolds 



A(t) := (1-tJAi+tAj 



(A39) 



as a function of < t < 1. Here, the initial, A i , and final, 
A f , configurations belong to G/H = U(2N r )/0(2N r ), 
while A(t) is not a member of G/H = U(2N r )/0(2N r ) for 
< t < 1. According to Eq. (A2S), the spectrum X L (t) 
of D[A(t)] is symmetric about the band center at the en- 
ergy zero. Configurations A. t and A^ have Pfaffians of 
opposite signs whenever an odd number of level crossing 
occurs at the band center ( "spectral flow" ) during the t- 
cvolution of the kernel D[A(t)]. This is accompanied by 
the closing of the spectral gap of D[A(t)] by an odd num- 
ber of pairs ( — A, (t), +A t (t)) as t interpolates between 
and 1. The spectral i-evolution is obtained numerically 
using the regularization of the kernel D[A(t)] by choos- 
ing the family on the torus T 2 . In this way, the index i 
takes discrete values. In Fig. we show the evolution 
of the eigenvalues for A(t) interpolating between A /=0 
and A i=1 . Observe that in A ; the part responsible for 
the winding configuration Ri(9,<f>) is entirely localized 
in the sector of the first replica. Thus, when comput- 
ing the spectral flow, we can focus on this sector alone. 
Since level crossing at the band center takes place for a 
single pair of levels, we conclude that Pf(D[A ;=0 ]) and 
Pf(U[A ; _ 1 ]) di ffer by their sign. This supports numeri- 
cally Eq. ^3§. 



There are 10 target spaces G/H for the NLcrMs of rel- 
evance ia the 10 symmetry classes of Anderson localiza- 
tion. EStIj They encode 10 distinct patterns of symmetry 
breaking. These patterns have been exhaustively classi- 
fied within a supcrsymmetric approach by Zirnbauer in 
Ref. [l2| Each target superspace G/H is a Riemannian 
symmetric supermanifold that can be parametrized in its 
bosonic sector by the Riemannian symmetric manifold 



M B = M BB x M FF . 



(Bl) 



Here, M B is the direct product between a non-compact 
Riemannian symmetric manifold M BB originating from 
the boson-boson sector of the Riemannian symmetric su- 
permanifold and a compact Riemannian symmetric man- 
ifold M FF originating from the fcrmion-fermion sector of 
the Riemannian symmetric supermanifold. The target 
superspaces G/H and Riemannian symmetric manifolds 
M B = M BB x M FF relevant to this paper are: 

• G/H = GL(n|n) x GL(n|n) / 'GL(n|n) = GL(n|n) 
with the Riemannian symmetric manifolds 



M 
M 



BB = GL(n,C)/U(n), 



ff -= U(n), 
for the chiral symmetry class AIII, 



(B2) 
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G/H = GL(2n|2n)/OSp(2n|2n) with the Rieman- 
nian symmetric manifolds 



M BB = GL(2n,K)/0(2n) 
M FF = U(2n)/Sp(2n), 
for the chiral symmetry class BDI, 



(B3) 



G/H = GL(2n|2n)/OSp(2n|2n) with the Rieman- 
nian symmetric manifolds 



M 
M 



BB 



(B4) 



FF 



U*(2n)/Sp(2n), 

U(2n)/0(2n), 

for the chiral symmetry class CII, 

G/H = OSp(2n|2n)/GL(n|n) with the Riemannian 
symmetric manifolds 



7\/ BB = Sp(4n,R)/U(n), 
M FF = 0(2n)/U(n), 



(B5) 



for the BdG symmetry class D, 

G/H = OSp(4n|4n)/OSp(2n|2n) x OSp(2n|2n) 
with the Riemannian symmetric manifolds 



(B6) 



Sp(2n,2n)/Sp(2n) x Sp(2n), 
l ff = SO(4n)/SO(2n) x S0(2n), 
for the symplectic symmetry class AIL 

The compatibility of these target superspaces with 
the addition of a topological term in the corresponding 
NLerM is solely determined by the compact Riemannian 
symmetric manifold M FF : A topological term requires 
a non-trivial second homotopy group of M FF (e.g., CII 
and All). In this context, observe that the Riemannian 
symmetric supermanifolds for symmetry classes BDI and 
CII merely differ by the exchange of the boson-boson and 
fcrmion-fcrmion stabilizers 0(2n) and Sp(2n) [in that re- 
gard, it is convenient to view U*(2n) as a non-compact 
real subgroup of GL(2n,C)]. This small difference is of 
great consequences since the NLcM for the symmetric 
class BDI cannot be augmented by a topological term. 

For simplicity, we cons ide r symmetry class CII with 
n = l. According to Eq. (B4), the target superspace is 

GL(2|2)/OS P (2|2) « GL(2|2)/SL(1|2) (B7) 

whereby we used the isomorphism OSp(2|2) w SL(1|2). 
The projected CII target superspace obtained by quoti- 
cnting out the two diagonal generators of GL(2|2) is 

PSL(2|2)/OS P (2|2) rj PSL(2|2)/SL(1|2). (B8a) 



One must carry this projection on the non-c omp act and 
compact Riemannian symmetric manifolds (B4). This 
is done by quotienting out their R + and U(l) factors, 
respectively. Hence, the projected boson-boson Rieman- 
nian symmetric manifold is 



SU*(2)/Sp(2) w SU*(2)/SU(2) 



while the projected fermion-fermion Riemannian sym- 
metric manifold is 



SU(2)/0(2) ~ S 2 



(B8c) 



Appendix C: The supergroup PSL(n|n) 

Let n £ N be an integer and denote with smat(n|n) 
the set of all real (n\n) supermatrices M, i.e., matrices 
of the form 



M 



M BB M BF 

FB M FF 



(CI) 



where M BB and M FF are n x n reaZ-valued matrices 
while M BF and M FB are n x n Grassmann-valued ma- 
trices. We shall denote with I and J the (n|n) diagonal 
supermatrices 



I := diag (l ■ • ■ 1 1 • • • l) , 
J:=diag(l ••• 1 -1 •■■ -1) 



(C2) 



respectively. For any element M £ smat(n|n), the super- 
trace is defined by 



str M := tr M BB — tr M FF . 



Observe that 

trM 



3 - - (str MJ + str MI) 



(C3) 



(C4) 



tr M BF = - (str MJ~ str MI) , 



i.e., demanding that trM BB and trM BF both vanish is 
equivalent to demanding that strMJ and str MI both 
vanish. For any element M £ smat(n|n) with det M FF ^ 
or det M BB =^ the superdeterminant is defined by 

det (M RR - M RF M~ F M FR ) , 
sdetM:= v BB , — BF FF FBJ , (C5) 



detM 



FF 



det (M FF - M FF ,M R lM F!F ) , x 

sdet M := BB BF \ C6 

det M BB 

respectively. 

An obvious generalization of smat(n|n) is achieved 
through the complexification 

M -> M + iM', M, M' £ smat(n|n). (C7) 

Another one follows from the substitution 



smat(njn) — > smat(m|n) 



(C8) 



where th e su permatrices from the set smat(m|n) are of 
the form (CI) with the entries of the mx m matrix M BB 



and the n x n matrix A'I FF commuting numbers while 
the entries of the the m x n matrix M BF and the n x m 



(B8b) matrix M FB are anticommuting numbers. 



The following definitions apply to both real and com- 
plex supermatrices. The supergroup PSL(n|n) is con- 
structed from the supergroup GL(n|ra) C smat(njn) as 
follows. The supergroup GL(n|n) consists of all (n\n) 
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supermatrices for which both M BB and M FF are non- 
singular (i.e., have nonvanishing determinants) and with 
the matrix multiplication as the group operation. The 
supergroup GL(n|n) is not semisimple. It possesses the 
matrix sub-supergroup SL(n|n) that follows from restrict- 
ing the superdeterminants in GL(n|n) to one. The super- 
group SL(n|n) is also not semisimple, for it contains the 
(n\n) unit supermatrix I that commutes with all (n|n) 
supermatrices. All elements of SL(n|n) are generated 
through exponentiation of elements of the Lie superal- 
gebra sl(n|n) whereby any element of sl(n|n) is a (n\n) 
supermatrix of the form 



A' 



^fb X ff 



(C9) 



with X BB and X FF nxn real- valued matrices while M BF 
and M FB are nxn Grassmann- valued matrices with the 



vanishing supcrtracc 



str A = 0. 



(CIO) 



The supergroup PSL(n|n) is defined to be the factor- 
group SL(n|n)/R + (R, the set of positive real numbers) 
by which any two elements in sl(n|n) that differ by a 
multiple of the unit element I generate upon exponentia- 
tion the very same element of PSL(n|n). The supergroup 
PSL(n|n) is semisimple, an element of PSL(n|n) cannot 
be written as a supermatrix. 



Appendix D: The Lie group U*(2) 

Let n £ N. The Lie group U*(2n) is the set of matrices 
in GL(2n,C) that commutes with the linear transforma- 
tion 



tp : C 2 



\ z 2n j 



u + 



'an 



(Dl) 



\-z n ) 



where complex conjugation is denoted by *. It follows 
that the Lie algebra u*(2n) is the set of matrices in 
GL(2n,C) of the form 



Z l Z 2 

-zi z\ 



(D2) 



where Z x and Z 2 are any complex- valued nxn matrices. 

We now specialize to the Lie group U*(2) with the Lie 
algebra u*(2). Let A Gu*(2). There exist the complex 
numbers z 1 and z 2 such that 



A = 



(D3) 



Here, we have introduced the unit 2x2 matrix cr and 
the three Pauli matrices ((T 1 , <t 2 , cr 3 ). Evidently, U*(2) 
and GL(l,R)xSU(2) share the same Lie algebra, 

u*(2) «R®su(2). (D4) 

Locally, one thus have the isomorphism 

U*(2) «R + ® SU(2) (D5) 

where R + is the set of positive real numbers. 



Appendix E: NLcrM on homogeneous target spaces 
versus NLcrM on symmetric target spaces 

In this Appendix we briefly summarize a number of 
facts about NLcrM on target spaces which are homo- 
geneous spaces G/H. (For more detai ls w e refer the 
reader, for example, to the appendix of 106, and refer- 
ences therein.) 

The essential properties of a NLcrM whose target space 
is of the form of a coset space G/H, where G is a Lie 
group and H a Lie-subgroup, are as follows. 

If H is a maximal subgroup of G. then the NLcrM has 
precisely one coupling constant (the target space G/H is 
then what is known as a 'symmetric space'). If, on the 
other hand, there exists precisely one intervening sub- 
group H', i.e., 



H <zH' c G, 



(El) 



: Re z 1 cr + *I m z 2 a l + *-R- c z 2 °2 + z l a 3- 



then the NLcrM with target space G/H turns out to have 
precisely two independent coupling constants. If, in the 
latter case, we were to run the RG into the infrared (i.e., 
to large length scales), then one of the two coupling con- 
stants will disappear. In this case one ends up, asymptot- 
ically at long scales, with a NLcrM with one coupling con- 
stant, whose target space will be either G/H' or H'/H 
(both are, by assumption, symmetric spaces in the sense 
of maximal subgroups). The number of coupling con- 
stants for a NLcrM on a homogeneous space with more 
intervening subgroups increases accordingly. 

The relevance of the symmetric space NLcrM for the 
universal physics appearing at large length scales is thus 
simply a consequence of the RG. If we begin with a 
NLcrM on a general homogeneous space, the RG will se- 
lect, asymptotically at long length scales, a target space 
which is a symmetric space G/H where H is a max- 
imal subgroup of G. So, all NLcrM on coset spaces be- 
come NLcrM on symmetric spaces, asymptotically at long 
length scales. It is for this reason that they are the 
"stable" large-scale limits, and appear naturally, with- 
out fine-tuning. 

The appearance of a NLcrM (on a, in general, possi- 
bly homogeneous and not necessarily symmetric space) 
arises in a physical context from the general principle of 
symmetry breaking. The group G is the global symme- 
try group of the problem. The subgroup H (not neces- 
sarily maximal) characterizes the symmetries which are 
preserved when the global symmetry G is broken. 



2G 



In the situation discussed in this article, the global 
symmetry group is 

G = GL(2to|2to), m=l,2,---. (E2) 

When a global density of states (DOS) is generated in 
the Dirac fcrmion formulation of the theory, the resulting 
expectation value breaks the symmetry G but preserves 
the symmetry H C G. It follows from Ref. [35]that on the 
"dotted line" in our global phase diagram in Fig. || this 
subgroup is W = OSp(2m|2?7i). The manifold G/H' = 
GL(2m|2m)/OSp(2m|2m) is the symmetric space corre- 



sponding to symmetry class CII. On the boundary of the 
CII region the symmetry class is All, corresponding to 
a target manifold H'/H = OSp(2m|2m)/[OSp(m|m) x 
OSp(m|m)]. Thus, the entire lower half of our global 
phase diagram in Fig. |] can be described by the NLcM 
on the homogeneous space 

G/H = GL(2m|2m)/[OSp(m|m) x OSp(m|m)]. (E3) 

which has two coupling constants (besides the coupling 
of th e Gade term) . This is analogous to the discussion in 
106. 
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(5.2) from Ref. M should be multiplied by 27r. (The same 
beta function is denoted by f3 g in the present paper.) 
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Low Dimensional Applications of Quantum field Theory, 
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Study Series BiPhysics (Plenum Press, New York, 1997), 
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The stability of the vertical dotted line in region BDI 
of Fig. ^ follows from the analytical continuation g A — » 
~5a'. 9v -> -ftmmi and g M +tecm of these flows, 
as pointed out in Ref. p4. Similarly, the stability of the 
vertical dotted line in region CII of Fig. ^ follows from the 
analytical continuation g A — } —g a i, gv +<?i m ?n'i an d 

9m -> -teem' of these flows- 
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aries g% = g 2 - or the CII boundaries g'+ — g'l. The pro- 
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in Fig. ^(b). Given these finite slopes of the flows out of 
the boundaries of regions BDI and CII, the reader might 
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